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ABSTRACT 


This  thesis  examines  spectrometric  oil  analysis  data 
from  two  sources  in  an  attempt  to  formulate  a  statistical 
model  which  will  be  useful  in  monitoring  aircraft  engines 
in  the  Naval  Oil  Analysis  Program.   Initially,  experimental 
data,  gathered  for  an  Air  Force  study,  is  used  to  determine 
if  the  measurement  error  inherent  in  the  monitoring 
procedure  is  normally  distributed  and  if  correlations  exist 
between  measurements  for  different  wear  metals.   Based  on 
the  results  of  this  investigation,  a  study  is  made  of 
operational  data  from  Wright  reciprocating  engines  of  the 
R1820-82  model  type.   This  investigation  leads  to  the 
conclusion  that  a  multivariate  regression  model  is  useful 
in  estimating  the  parameters  of  the  distribution  of  analyses 
from  properly  operating  engines  of  this  type.   A  procedure 
is  then  suggested  which  would  employ  the  readings  from  past 
oil  analyses  from  a  particular  engine  to  determine  its 
present  condition. 
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I.   INTRODUCTION 

For  several  years  the  technique  of  using  the  spectro- 
metric  analysis  of  oil  samples  as  an  aid  in  determining  the 
condition  of  diesel  engines  has  been  employed  successfully 
by  major  railroads  and  various  other  users  of  large  diesel 
equipment.   In  1956,  a  trial  program  was  begun  at  the  Naval 
Air  Rework  Facility  in  Pensacola  to  determine  if  this 
method  could  also  be  used  to  monitor  aircraft  engines. 
Since  that  time  the  program  has  proved  successful  and  has 
evolved  into  the  Naval  Oil  Analysis  Program  (NOAP) .   It  is 
planned  that  this  program  will  eventually  include  all  Navy 
fluid  lubricated  mechanical  systems.   A  more  detailed 
history  of  the  program  is  contained  in  Refs.  1  and  2. 

Since  this  thesis  is  concerned  with  an  investigation  of 
data  collected  at  the  Pensacola  laboratory,  the  following 
descriptions  of  the  operation  will  be  limited  to  the  pro- 
cedures used  there.   Reciprocating  aircraft  engines  are 
sampled  approximately  every  30  hours.   The  sample  is  taken 
after  the  aircraft  has  returned  from  a  flight  and  before 
the  oil  has  become  cold.   It  is  immediately  sent  to  the 
laboratory  by  air  mail  and  is  analyzed  on  the  day  received. 
The  analysis  is  accomplished  by  a  spectrometer  using  the 
rotating  graphite  electrode  technique.   Measurements  of  the 
parts  per  million  (ppm)  content  of  ten  metallic  elements, 
which  might  be  indicative  of  engine  wear,  are  made 


simultaneously.  Of  these  ten,  aluminum,  chromium,  iron, 
silver,  copper,  magnesium  and  nickel  are  those  which  are 
relevent  to  engines  of  the  model  considered  in  this  report. 
The  ppm  readings  are  automatically  recorded  on  a  punched 
card  which  also  contains  various  other  hand-entered  data 
identifying  the  sample. 

Once  the  data  has  been  recorded,  it  is  used  to  aid  in 
determining  what  the  operating  condition  of  the  engine 
might  be.   Presumably,  if  the  engine  is  operating  properly, 
the  amount  of  metallic  contamination  in  the  circulating  oil 
should  be  within  certain  normal  limits.   In  addition,  it  is 
felt  that  the  amount  of  contamination  added  to  the  oil 
since  it  was  last  sampled  should  be  within  specific  limits 
if  the  engine  is  in  good  working  condition.   If,  however, 
the  engine  is  discrepant  and  excessive  wear  is  present, 
this  will  presumably  cause  an  abnormal  addition  of  metallic 
contaminates  to  the  circulating  oil. 

Thus,  when  a  sample  has  been  analyzed  and  the  results 
recorded,  both  the  magnitudes  of  the  present  readings  and 
the  changes  in  the  readings  since  the  last  sampling  are 
compared  with  threshold  limits  which  have  been  developed 
for  each  engine  type  and  each  metallic  element  relevent  to 
that  engine.   If  the  results  fall  outside  the  prescribed 
limits,  some  action  is  generally  taken  by  the  laboratory. 
Usually  another  sample  is  requested  and  the  previous 
results  are  verified.   If  the  abnormality  persists,  either 


the  aircraft  is  grounded  for  maintenance  or  future  samples 
are  taken  more  frequently  than  the  usual  30  hour  interval. 

At  present,  these  threshold  limits  are  subjectively  set 
and  vary  only  from  element  to  element  and  among  engine 
model  types.   They  are  based  on  the  past  history  of  the 
aircraft  model  which  includes  the  data  supplied  by  the 
engine  manufacturer  before  the  model  is  placed  in  service 
and  experience  accumulated  once  the  model  is  in  use.   The 
limits  are  not  used  as  sharp  boundaries  for  classifying 
engines  as  normal  or  discrepant  but  merely  as  indicators 
upon  which  a  subjective  decision  as  to  the  action  to  be 
taken  can  be  based. 

This  report  examines  two  sets  of  data  from  the  Pensa- 
cola  laboratory  with  the  intention  of  determining  the 
propriety  of  three  assumptions  implicit  in  this  classifi- 
cation procedure.   Since  the  same  threshold  limits  are 
used  for  all  engines  of  a  particular  model,  it  is  assumed 
that  all  normally  operating  engines  of  the  same  type  can  be 
expected  to  have  the  same  amounts  of  metallic  contamination 
in  their  oil  systems.   In  addition,  since  threshold  limits 
are  constant  for  a  given  element  and  model  type,  variations 
in  other  factors,  such  as  the  operating  hours  since  the 
last  oil  change,  must  be  ignored  or  subjectively  introduced 
into  the  classification  procedure.   Finally,  since  thresh- 
old limits  are  set  for  each  element  independent  of  the 
limits  for  other  elements,  readings  for  different  elements 
are  assumed  to  be  uncorrelated. 


Once  these  assumptions  are  verified  or  rejected,  a 
statistical  model  is  formulated  to  aid  in  establishing  a 
more  objective  classification  criterion. 


II.   ERRORS  INHERENT  IN  THE  MONITORING  PROCEDURE 

Since  the  intention  of  NOAP  is  to  make  inferences 
about  the  condition  of  aircraft  engines,  based  on  the 
amount  of  wear  metal  contamination  in  the  engine's  oil 
system,  it  is  extremely  important  that  the  amount  of  con- 
tamination recorded  at  the  laboratory  accurately  reflect 
the  actual  amount  present  in  the  engine.   For  the  purposes 
of  this  report,  measurement  error  will  be  defined  as  the 
difference  between  the  parts  per  million  content  of  an 
element  recorded  as  present  in  a  particular  engine  at  a 
point  in  time  and  the  actual  content  at  that  time.   In  NOAP 
there  are  a  variety  of  potential  sources  of  error,  all  of 
which  can  contribute  to  the  net  measurement  error  defined 
above . 

A.   ERRORS  IN  SAMPLING 

As  was  mentioned  earlier,  oil  samples  are  taken  from 
reciprocating  engines  normally  every  30  flying  hours  and 
while  the  oil  is  still  hot.   This  sampling  is  accomplished 
with  a  special  sampling  kit  consisting  of  a  sampling  tube 
and  a  sampling  bottle.   The  tube  is  inserted  into  the  oil 
reservoir,  and  when  it  has  filled  the  top  end  is  stopped 
with  the  operator's  finger.   The  contents  are  then  trans- 
ferred to  the  bottle,  which  is  immediately  forwarded  to 
the  laboratory  for  analysis.   When  the  sample  is  analyzed, 
a  small  portion  of  the  oil  in  the  bottle  is  used  in  the 


analysis  [Refs.  1  and  2].   Thus,  an  extremely  small  amount 
of  oil  is  used  to  determine  the  extent  of  contamination  in 
the  engine's  entire  oil  system.   Any  lack  of  homogeneity 
in  the  engine's  oil  reservoir  will  result  in  a  non- 
representative  sample.   Further,  any  contamination  added 
to  the  sample  through  a  lack  of  cleanliness  of  the  sampling 
tube  and  bottle  or  through  handling  at  the  laboratory  will 
contribute  to  the  measurement  error. 

B.  ERRORS  IN  RECORDING 

At  the  time  the  sample  is  taken,  certain  data  including 
the  date,  the  operating  hours  since  the  last  oil  change, 
the  hours  since  the  last  overhaul  of  the  engine,  the  engine 
serial  number  and  the  model  number  are  recorded  and  mailed 
to  the  laboratory  with  the  sample.   Various  portions  of 
this  data  are  transferred  from  other  records.   At  the 
laboratory  the  data  are  entered  by  hand  on  the  permanent 
record  cards  maintained  there  [Ref.  1] .   This  entire 
sequence  of  recording  and  transferring  data  from  one  record 
to  another  can  result  in  mistakes. 

Unfortunately,  as  with  the  errors  in  sampling,  there  is 
no  data  available  at  the  present  time  that  can  be  used  to 
measure  this  error. 

C.  ERRORS  IN  ANALYSIS 

When  an  oil  sample  is  received  at  the  Pensacola  labora- 
tory, it  is  analyzed  using  a  direct  reading  spectrometer 
with  spark  excitation,  stationary  and  rotating  disc 
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electrodes.   The  sample  bottle  cap  is  filled  with  oil  and 
placed  in  the  spark  stand.   The  gap  between  the  two 
electrodes  is  set  and  the  disc  electrode  begins  to  rotate 
at  30  rpm.   As  the  electrode  rotates,  a  thin  film  of  oil 
is  forced  to  the  area  under  the  fixed  electrode.   A  high 
energy  spark  is  then  fired  across  the  gap  and  the  film  of 
oil  is  burned  for  25  seconds.   The  light  from  the  burning 
oil  is  separated  so  that  its  intensity  at  the  wave  lengths, 
produced  by  the  elements  to  be  analyzed,  can  be  compared 
with  built-in  standards.   The  average  intensity  over  the 
burning  period  is  then  measured  for  each  element  simul- 
taneously and  converted  into  parts  per  million.   These 
readings  are  automatically  recorded  on  the  engine  history 
card  [Ref .  2] . 

If  it  is  assumed  that  there  were  no  errors  in  the 
sampling  or  recording  and  thus,  that  the  oil  used  in  the 
analysis  is  representative  of  the  oil  in  the  engine,  any 
difference  between  the  true  content  of  contamination  in 
the  engine  and  that  recorded  after  the  analysis  can  be 
attributed  to  an  analysis  error.   An  experiment  designed 
to  measure  the  effects  of  this  type  of  error  has  been 
conducted  and  the  results  of  an  analysis  of  this  data  are 
presented  in  the  next  section. 
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III.   DISCUSSION  OF  THE  ANALYSIS  ERROR 

Although  the  error  due  to  the  spectrometric  analysis 
of  the  oil  is  not  the  only  possible  source  of  error,  it 
certainly  is  a  major  contributor  to  the  over-all  measure- 
ment error  defined  earlier.   For  this  reason,  an  examination 
of  data  accumulated  for  a  study  conducted  by  the  Air  Force 
[Ref.  3]  was  performed  and  the  results  are  discussed  in 
this  section. 

A .   DATA 

In  1967,  the  Pensacola  laboratory  participated  in  an 
experiment  conducted  by  the  Air  Force.   During  a  30  day 
period  the  laboratory  received  100  oil  samples.   These 
were  to  be  analyzed  in  the  normal  manner  and  the  results 
reported.   Although  the  laboratory  was  not  aware  of  it, 
these  100  samples  consisted  of  ten  samples  each  repeated 
ten  times.   Thus,  the  laboratory  actually  repeated  the 
analysis  of  ten  different  samples  ten  times.   The  results 
of  the  analyses,  for  the  seven  elements  of  interest  in 
this  report,  are  included  in  the  Computer  Output  section, 
where  the  readings  on  like  samples  are  in  groups  numbered 
from  one  to  ten.   Missing  data  accounts  for  some  groups 
having  less  than  ten  readings. 

For  each  of  the  ten  groups  of  repetitious  analyses,  the 
sample  mean  and  standard  deviation  were  calculated  for  each 
element.   If,  for  example,  X.  is  the  i —  reading  for 
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aluminum  in  a  group  of  size  n,  then  the  sample  mean,  X, 
and  standard  deviation,  S,  for  aluminum  are 


n 
X  =  (  £  X. )/n 
i-1  X 


and 


n       _  2       h 
S  =  [  E  (X.  -  X)  /(n-1)]  (1) 

i=l   x 


respectively.  The  results  of  these  calculations  are  also 
presented  in  the  Computer  Output  section  for  each  element 
and  each  group  of  samples. 

In  addition,  for  each  sample  group,  an  estimated  corre- 
lation matrix,  R,  was  calculated,  where  if  r. .  is  the 

—  lj 

element  in  the  i —  row  and  j —  column  of  R,  then 

n 

Z  (X.  .  -  X. ) (X.  .  -  X.) 

r.  .  =  £=± j-        (2) 

ID      n  _   n  9  h 

[    I     (X.  ,  -  X. )Z  Z     (X.  .  -  X.)Z] 
k=l   ^k     X    k=l   ^'k     1 

where  X.  ,  is  the  k —  reading  for  the  i —  element,  and  X. 

1  ,  K  1 

is  the  element's  sample  mean.  These  correlation  matrices 
are  included  in  the  Computer  Output  section. 

These  preliminary  computations  provided  statistics 
which  were  used  to  test  certain  hypotheses  concerning  the 
probability  distribution  of  the  analysis  error. 
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B.   TEST  FOR  NORMALITY 

Since  the  overall  measurement  error  is  the  net  effect 
of  errors  arising  from  a  variety  of  sources,  the  Central 
Limit  Theorem  of  Probability  Theory  [Ref.  4]  provides  good 
justification  for  making  an  assumption  of  normality  in  the 
distribution  of  this  error.   Thus,  any  additional  evidence, 
which  tends  to  point  to  the  normality  of  one  of  the  con- 
tributing sources  of  error,  will  serve  to  strengthen  the 
overall  assumption.   The  experiment  conducted  by  the  Air 
Force  provided  data  which  was  used  to  test  for  normality 
in  the  distribution  of  the  analysis  error.  ■  •• 

If  X,  is  defined  as  a  seven-component  vector  of  sample 
readings  arising  from  sample  group  k  =  1,2,...,  10,  then 


^k  =  Hk  +  £* 


where  y,  is  a  seven-component  vector  of  the  true  metallic 

content  of  seven  elements  in  the  sample  associated  with  the 

k —  group  of  readings,  and  e,  is  the  seven-component  random 

analysis  error  vector.   Thus,  if  it  is  assumed  that  e,  is 

a  multivariate  normal  random  variable  with  zero  mean  vector 

and  unknown  covariance  matrix  Z,  ,  then  X,  is  a  multivariate 

normal  random  variable  with  mean  vector  y,  and  covariance 

— k 

matrix  Z_,  ,  denoted  N  (y.  ,  £_,  )  . 

If  this  assumption  is  correct,  it  is  possible  to  perform 
a  transformation  of  the  form, 


Z.  =  P.  (X,  -  y,  )  , 

— k    — k  — k    — k 
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which  will  produce  a  multivariate  normal  random  variable  Z, 
which  has  mean  vector  zero  and  covariance  matrix  I,  the 
identity  matrix.   For  details  of  this  transformation  see 
Appendix  A. 

For  each  of  the  ten  groups  of  sample  readings  the  mean 
vector  |jl  was  estimated  using 


nk 


4  =  .z  xkj/"k 

1  =  1   J 


(3) 


where  X,  .  is  the  j —  vector  of  sample  readings  from  the  k — 
sample  group,  and  n,  is  the  number  of  readings  in  that 
group.   In  addition,  the  covariance  matrices,  E_,  ,  were 
estimated  using  the  unbiased  estimator 


E,  =  (S.  S  .  r.  .),  (4) 

%     i   3   ij  k 

where  S.  is  the  estimated  standard  deviation  for  the  i — 

i 

element,  computed  as  in  equation  (1) ,  and  r. .  is  as  defined 
by  equation  (2) .   For  each  of  the  ten  groups  the  non- 
singular  matrix,  P ,  ,  was  found  and  the  transformation, 

Z,   .  =  P.  (X,   .  -  X.  )  , 
— k  ,  l    — k  — k ,  l    — k  ' 

4- Vi 

performed  on  each  vector  of  readings,  X,  . ,  in  the  k — 
group,  k  =  1 , 2 , . . . ,10 .   In  this  way  vectors  Z,  .  were 

— K  ,  1 

produced,  the  components  of  which  are  stochastically  inde- 
pendent and  are  distributed  according  to  N (0,1)  if  the 
hypothesis  is  true. 
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All  readings  from  the  ten  groups  were  then  pooled  to 
produce  a  sample  of  651  deviates,  assumed  to  be  univariate 
normal.   The  Kolmogorov-Smirnov  goodness  of  fit  test  was 
applied  to  this  sample  and  the  resulting  test  statistic  of 
.0215  was  not  significant  even  at  the  .20-level.   Thus,  the 
hypothesis  of  normality  in  the  distribution  of  the  analysis 
error  was  accepted.   For  details  of  the  test  see  Appendix  A, 

C.   TEST  FOR  EQUALITY  OF  COVARIANCE  MATRICES 
Let  X,  again  be  defined  as 

*k  =  Hk  +  ^k 

as  in  the  previous  section,  where  now  e,  is  assumed  to  be 
a  multivariate  normal  random  vector.   In  addition,  let  the 
estimators  X,  and  E,  be  defined  by  equations  (3)  and  (4) 
respectively.   Then,  the  Air  Force  data  can  be  used  to 
produce  ten  estimated  covariance  matrices  E,  ,  each  asso- 
ciated with  a  different  sample  group  and  thus,  a  different 

true  content  vector  un  .   If  the  true  covariance  matrices, 

— k 

E,  ,  are  independent  of  the  vector  y_ ,  ,  and  thus  constant 
for  all  k  =  1,2 ,...,10,  the  ten  estimated  matrices  could 
be  pooled  to  obtain  an  over-all  estimate  of  E.   This 
hypothesis  of  the  equality  of  the  ten  covariance  matrices 
was  tested,  and  the  results  led  to  the  rejection  of  the 
hypothesis  at  the  .10  level  of  significance.   The  details 
of  the  test  used  and  the  results  obtained  are  included  in 
Appendix  A. 
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In  an  attempt  to  account  for  the  apparent  variability 
among  covariance  matrices  from  different  sample  groups,  a 
regression  model  was  formulated.   It  had  been  suggested  by 
Baird-Atomic  Inc. ,  the  manufacturer  of  the  spectrometer 
used  at  the  Pensacola  laboratory,  that  the  variability  in 
readings  for  a  given  element  is  dependent  on  the  true 
content  of  the  element.   Specifically,  the  relationship 
is  assumed  to  be 

2       ,2 

a   =  a  +  by 

2  . 

where  a   is  the  variance  of  repeated  analyses  of  the  same 

sample  for  a  given  element,  y  is  the  true  parts  per  million 
content  of  the  element,  and  a  and  b  are  constants  [Ref.  5]. 
The  model, 

S2  =  a  +  bX2  +  e, 

was  used  to  examine  the  propriety  of  this  relationship  for 

2 
each  of  the  seven  elements  under  consideration.   Here  S 

is  the  square  of  the  standard  deviation  estimate  defined  by 

_2 
equation  (1),  X   is  the  square  of  the  sample  mean,  and  e  is 

a  random  variable.   For  each  element,  a  and  b  were  estimated 

using  least-squares  techniques  and,  under  the  assumption 

that  variations  about  the  regression  line  are  normally 

distributed,  the  hypothesis,  b  =  0,  was  tested  [Refs.  6 

and  7] .   A  t-test  was  used  and  the  significance  level  set 

at  .10.   Of  the  seven  slopes  tested,  those  associated  with 

aluminum,  iron,  copper  and  magnesium  were  significantly 
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non-zero.   It  should  be  mentioned  that  the  true  content  of 
the  other  elements  did  not  vary  much  among  the  ten  sample 
groups.   The  results  of  the  least-squares  estimation, 
together  with  the  numerical  results  of  the  t- tests,  are 
included  in  the  Computer  Output  section. 

D.   TEST  FOR  INDEPENDENCE  AMONG  ELEMENTS 

Because  of  the  apparent  dependence  of  the  variance  of 
repeated  readings  upon  the  true  content  of  an  element,  it 
was  felt  that  covariances  between  elements  might  depend 
upon  the  content  of  the  elements  concerned.   If  this  were 
the  case,  then,  for  example,  a  particularly  high  reading 

- 

of  one  element  might  be  "explained"  by  a  corresponding 
low  reading  of  another  element.   Under  the  assumption  of 
normality,  the  hypothesis  of  independence  is  equivalent  to 
the  hypothesis  of  zero  correlation.   This  hypothesis  was 
tested,  using  each  of  the  ten  correlation  matrices  R, 
defined  by  equation  (2) .   Of  the  ten  tests  conducted,  only 
three  of  them  were  not  significant  at  the  .10  level.   The 
numerical  results  and  the  test  used  are  given  in  Appendix  A. 
It  can  be  noted  from  the  correlation  matrices  in  the 
Computer  Output  section  that  particularly  strong  corre- 
lations seem  to  exist  between  iron  and  copper,  silver  and 
copper,  and  magnesium  and  iron. 

Thus,  it  appears  that,  in  general,  the  readings  of 
different  elements  are  not  independent,  and  some  explana- 
tion, for  example,  of  an  erroneous  copper  reading  may  come 
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from  an  examination  of  the  corresponding  iron  reading  on 
the  same  sample. 
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IV.   ANALYSIS  OF  OPERATIONAL  DATA 

Using  the  evidence  provided  by  the  Air  Force  data  to 
support  the  assumption  of  normality  in  the  distribution 
of  the  measurement  error,  an  investigation  of  some  opera- 
tional data  was  conducted.   A  statistical  model,  which  makes 
use  of  the  apparent  correlations  between  the  readings  on 
different  elements,  was  formulated.   The  details  of  the 
model's  formulation  are  discussed  in  this  section. 

A.   NOAP  DATA 

The  Naval  Air  Rework  Facility  at  Pensacola  provided  a 
magnetic  data  tape  containing  the  records  of  operational 
analyses  performed  there  from  the  beginning  of  July  to  the 
end  of  September  in  1967.   The  records  of  some  21,000 
different  analyses  were  included  on  the  tape.   For  each 
analysis  the  engine  model  number  and  serial  number,  as  well 
as  the  date  the  analysis  was  performed  and  its  results  in 
parts  per  million  for  each  relevant  element,  are  recorded. 
In  addition,  it  includes  the  operating  hours  since  the  last 
oil  change  and  since  the  last  overhaul  of  the  engine  for 
each  sample.   Unfortunately,  the  action  recommended  by  the 
laboratory  after  each  sample  was  analyzed  and  the  results 
of  that  action  were  not  available  with  the  tape.   For  this 
reason,  there  was  no  way  of  determining  with  certainty 
which  analyses  were  on  oil  from  properly  operating  engines 
and  which  were  from  discrepant  engines. 
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Of  the  113  different  engine  models  represented  on  the 
tape,  the  Wright  reciprocating  engine  model,  R1820-82,  was 
selected  for  investigation  since  it  was  the  most  frequently 
sampled,  with  4,134  different  analyses. 

B.   PRELIMINARY  RESULTS 

Since  it  seemed  logical  to  expect  the  content  of 
metallic  contamination  to  show  some  increase  from  normal 
wear  in  a  properly  operating  engine  as  the  hours  since  the 
last  oil  change  increase,  some  preliminary  plots  were  made 
by  the  computer.   Six  hundred  different  analyses  were  used 
with  no  regard  to  the  particular  engine  of  the  R1820-82 
type  from  which  they  came.   For  each  of  seven  elements, 
relevant  to  the  monitoring  of  engines  of  this  type,  the 
computer  plotted  the  ppm  content  versus  the  operating  hours 
since  the  last  oil  change.   For  three  of  the  elements, 
iron,  copper  and  aluminum,  there  was  some  indication  that  a 
buildup  of  contamination  occurs.   The  other  four  plots  gave 
no  evidence  of  any  significant  trend.   For  comparison,  the 
plot  of  chromium  is  included  with  those  of  aluminum,  iron, 
and  copper  in  Figures  1  to  4 ,  respectively.   Since  these 
plots  were  made  on  sample  readings  from  a  variety  of  en- 
gines, they  did  not  indicate  whether  a  particular  engine 
can  be  expected  to  show  the  same  trend.   For  this  reason, 
the  five  most  frequently  sampled  engines  of  the  R1820-82 
type  were  selected  and  the  computer  was  again  used  to  plot 
the  data  from  these  engines.   Five  different  symbols, 
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x,  +,  A,  Q  ,    and  D  ,  were  used  on  each  plot  to  represent 
the  five  different  engines.   Thus,  both  the  behavior  of  the 
readings  for  a  given  engine  and  differences  among  the 
readings  of  the  engines  could  be  seen.   Again,  evidence 
of  trends  was  limited  to  the  three  elements  iron,  copper 
and  aluminum.   For  these  elements  each  of  the  five  engines 
showed  a  roughly  linear  increase  in  the  ppm  content  as  the 
hours  since  the  last  oil  change  increased.   For  comparison, 
the  plot  for  chromium  is  again  included  with  the  plots  of 
aluminum,  iron  and  copper  in  Figures  5  through  8. 

C.   REGRESSION  MODEL 

Because  of  the  evidence  of  a  linear  increase  in  ppm 
content  versus  an  increase  in  hours  since  oil  change, 
provided  by  the  computer  plots,  a  regression  model  was 
suggested.   With  this  type  of  model  the  expected  content  of 
a  metallic  element  would  change  as  the  hours  since  the  last 
oil  change  varies.   Thus,  differences  between  what  is  a 
normal  amount  of  contamination  for  a  properly  operating 
engine  just  after  its  oil  has  been  changed  and  several 
flying  hours  later  could  automatically  be  incorporated 
into  a  classification  criteria. 

1.   Data  Selection 

Of  all  the  engines  of  the  R1820-82  model  type 
represented  on  the  tape,  those  with  eight  or  more  different 
analyses  were  selected.   Of  these,  any  with  missing  data  on 
one  or  more  records,  which  brought  the  usable  number  of 
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readings  to  less  than  eight,  was  rejected.   The  remaining 
engines  were  then  screened  in  an  attempt  to  insure  that  all 
the  data  had  come  from  properly  operating  engines.   As  was 
mentioned  earlier,  the  tape  did  not  provide  data  on  either 
the  recommendation  made  by  the  laboratory  or  the  results 
of  any  maintenance  which  might  have  been  recommended.   For 
this  reason,  it  seemed  that  the  best  way  to  insure  that  none 
of  the  data  had  come  from  a  discrepant  engine  was  to  exclude 
all  readings  which  exceeded  the  threshold  limits  presently 
used.   These  limits  are  given  in  Ref.  1.   After  this  pro- 
cedure had  been  applied  to  the  data,  any  engine  with  less 
than  eight  readings  was  eliminated  from  further  examination. 
At  this  point  there  were  27  engines  of  the  R1820-82  type 
remaining,  each  of  which  was  represented  with  from  eight 
to  fourteen  readings.   This  was  the  data  that  was  used  in 
the  remainder  of  the  investigation. 

2 .   Estimation  and  Tests  of  Regression  Coefficients 
It  was  assumed  that  the  outcome  of  each  spectro- 
metric  analysis  on  a  given  engine  is  of  the  form 

Y  =  BX  +  e 

where  Y  is  a  seven-component  vector  of  ppm  metallic 
contents;  B  is  a  7x2  matrix  of  unknown  coefficients;  X  is 
a  two-component  vector  with  first  component  identically 
one  and  second  component  equal  to  the  operating  hours  since 
the  last  oil  change;  and  e  is  a  multivariate  normal  error 
vector  with  mean  vector  zero  and  unknown  covariance  matrix, 
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E.   It  will  be  recalled  that  the  examination  of  the  Air 
Force  data  in  section  III  indicated  that  the  covariance 
matrix  may  vary  as  the  mean  content  vector  changes.   How- 
ever, the  slopes  of  the  regression  lines  mentioned  in  part 
C  of  that  section  are  such  that  variations  in  the  mean 
vectors  of  the  extent  present  in  the  operational  data  do 
not  result  in  an  appreciable  variation  in  the  elements  of 
E.   For  this  reason,  it  will  be  assumed  that  the  covariance 
matrix,  E_,  is  constant  in  the  development  to  follow.   Under 
this  assumption,  the  matrix  B  associated  with  each  engine 
was  estimated  using  the  least-squares  estimation  technique 
described  in  Appendix  B. 

Since  the  random  error  vector  is  assumed  to  be 
N  (0_,  z)   ,  each  observation  from  a  particular  engine  is 
N(BX,  E).   If  the  matrix  B  is  partitioned  so  that 


§.  =  (B^  B_2)  , 


then  a  test  of  the  value  of  each  of  the  components  of  B« 
can  be  made  to  determine  whether  the  variability  of  the 
readings  for  a  specific  element  is  related  to  variations 
in  the  hours  since  oil  change.   The  details  of  this  test 
are  included  in  Appendix  B.   For  each  engine  and  each 
component  of  B„  a  test  of  the  hypothesis  that  the  component 
is  equal  to  zero  was  made.   Table  1  gives  the  number  of 
times  a  specific  component  was  significantly  positive  or 
negative  at  an  over-all  a-level  of  .10.   For  this  a  value, 
the  expected  number  of  times  in  2  7  tests  the  results  will 
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be  significantly  positive  or  negative,  if  the  component  is 
actually  zero,  is  1.35. 

TABLE  1 

NUMBER  OF  SIGNIFICANT  REGRESSION  SLOPES 
TWO-TAILED  TEST;  a  =  . 10 


*2 


Element 


Number  of 
significantly 
positive 
components 


Number  of 
significantly 
negative 
components 


B 


B 


B 


B 


B 


B 


B 


2,1 
2,2 
2,3 
2,4 
2,5 
2,6 


2,7 


Aluminum 

Iron 

Chromium 

Silver 

Copper 

Magnesium 

Nickel 


8 
17 
3 
1 
16 
5 
5 


0 
0 
1 
3 
1 
0 
0 


These  results  indicate  that  the  metallic  content  of 
aluminum,  iron  and  copper  in  properly  operating  engines  of 
the  R1820-82  type  tends  to  increase  as  the  hours  since  the 
engine's  oil  was  last  changed  increase.   The  evidence 
pointing  to  this  conclusion  is  particularly  strong  in  the 
case  of  iron  and  copper.   Further,  there  seems  to  be  no 
significant  indication  that  such  a  relationship  exists 
in  general  for  chromium,  silver,  magnesium  or  nickel.   The 
numerical  results  of  the  tests  of  the  components  of  B_  for 
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each  engine,  in  addition  to  the  raw  data  and  estimates  of 
B  and  £_,  are  included  in  the  Computer  Output  section. 

Any  use  of  these  results  in  establishing  an 
operational  classification  procedure  for  identifying  dis- 
crepant engines  depends  upon  the  estimation  of  the  unknown 
matrix,  B,  from  past  analyses.   For  this  reason,  it  is 
important  to  determine  if  the  observations  from  different 
engines  all  come  from  the  same  over-all  probability 
distribution.   If  this  is  the  case,  all  data  on  engines  of 
the  R1820-82  model  type  could  be  used  to  estimate  a  single 
matrix,  B.   As  a  first  step  in  this  direction,  the  model 


Y.  =  B.X  +  e. 
—l    —l—   —l 

was  used  where  Y.  is  a  seven-component  vector  of  readings 


on 


engine  i;  B.  =  (Bn ,  B„) . ,  a  7x2  matrix  of  coefficients 
^        —l     —1   — 2  l 

associated  with  the  i —  engine  and  where  the  components  of 

B_  associated  with  chromium,  silver,  magnesium  and  nickel 

are  assumed  to  be  zero;  and  e.  is  N(0,  2.).   The  elements 

—l      —  —l 

of  B . ,  not  assumed  to  be  zero,  were  estimated  as  before 
—l 

and  used  to  estimate  the  2  7  covariance  matrices,  E..   The 

—l 

unbiased  estimate  of  £  .  is 

—l 

l.    =  ~   Z   (Y.    -  B.X.)  (Y.    -B.X.)' 

-l   ni-2  .  x  -i,:   -i-3  -i,:   -1-: 

th 
where  Y.    is  the  j —  observation  of  the  vector  Y . ;  B.  is 
-i,D         J  -i   -i 

j_t_ 

the  estimate  of  B.;  X.  is  the  j —  observation  of  X;  and 

— l   —  j  J  — 

n.  is  the  number  of  observations  associated  with  the  i — 
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engine.   A  test  of  the  hypothesis  of  equal  covariance 
matrices  was  made  using  these  estimates.   The  test  statistic 
was  extremely  significant  at  the  .10  level,  and  the  hypo- 
thesis of  equal  covariance  matrices  was  rejected.   The  test 
used  and  numerical  results  are  included  in  Appendix  A. 

Since  the  evidence  indicates  that  the  covariance 
matrices  associated  with  readings  from  different  engines 
are  not  the  same,  the  overall  conjecture  of  like  distri- 
butions must  also  be  rejected. 
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V.   CONCLUSION 

The  investigation  of  the  Air  Force  data  and  the  actual 
analysis  records  of  a  three  month  period  from  the  Pensacola 
laboratory  lead  to  three  main  conclusions.   First,  the  error 
inherent  in  the  ppm  spectrometric  readings  is  multivariate 
normally  distributed  with  significant  covariances  existing 
between  the  readings  of  various  pairs  of  elements.   Further, 
there  appears  to  be  a  linear  increase  in  the  content  of 
aluminum,  copper  and  iron  present  in  properly  operating 
engines  of  the  R1820-82  type  as  the  hours  since  the  last 
oil  change  increase.    Finally,  there  seems  to  be  no  justi- 
fication for  expecting  readings  on  samples  from  different 
engines  of  the  R1820-82  type  to  vary  in  the  same  manner. 
Based  on  these  results,  an  objective  classification 
criterion  can  be  formulated  which  may  be  of  use  in  imporv- 
ing  the  present  classification  procedure. 

For  example,  all  back  data  on  a  particular  engine  of 
the  R1820-82  type  which  was  in  proper  working  order  could 
be  used  to  estimate  the  matrix  B,  which  in  turn  could  be 
used  to  estimate  the  covariance  matrix  E_.   Then,  any 
observation,  Y,  of  the  spectrometric  analysis  of  a  new  oil 
sample  from  that  engine  is  distributed  as  N (BX,  L)    if  the 
engine  is  operating  properly.   The  estimates  of  B  and  !_ 
could  be  used  to  construct  a  confidence  region  R  (x) 
[Appendix  B] .   The  region  R  (x)  would  be  constructed  so 
that,  if  the  engine  is  operating  properly,  the  reading  Y 


36 


will  be  contained  in  the  region  R  (x)  with  probability  1-a. 
Thus,  one  classification  criterion  would  be:   classify  the 
engine  as  operating  properly  if  Y_  is  within  R  (x)  and 
classify  as  discrepant  otherwise.   By  making  a  small,  say 
.01,  the  number  of  operational  engines,  which  are  mistakenly 
classified  as  discrepant,  can  be  expected  to  be  of  the  order 
of  1  in  100.   However,  the  smaller  the  parameter  a  is  made, 
the  larger  the  region  R  (x)  becomes,  and  thus,  the  more 

LA 

likely  it  is  that  a  discrepant  engine  will  be  classified 

as  operating  properly. 

For  this  reason,  it  may  be  more  appropriate  to  use 

two  values  of  a.   For  example,  a  could  be  set  at  .10  and  a1 

at  .01  and  two  regions  R  (x)  and  R  , (x)  constructed.   In 

this  way  a  procedure  could  be  used  which  would  1)  classify 

the  engine  as  in  proper  working  order  if  Y  is  in  R  (x) ; 

2)  classify  as  discrepant  if  Y  is  not  in  R  , (x) ;  and  3) 

require  verification  of  Y  or  more  frequent  sampling  if  Y 

is  in  R  , (x)  but  not  in  R  (x) . 
a  a 

The  final  selection  of  a  specific  classification 
criterion  and  the  setting  of  the  appropriate  level (s)  of 
a  must  be  done  subjectively  and  should  be  based  upon  an 
examination  of  the  costs  involved.   If  the  cost  of  classi- 
fying a  discrepant  engine  as  operational  is  much  larger 
than  the  cost  of  grounding  an  operational  aircraft  then  a 
should  be  made  appropriately  large  compared  to  its  value 
if  the  reverse  were  true. 
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APPENDIX  A 

STATISTICAL  TESTS 

A.   TRANSFORMATION  OF  N (y ,  Z)  TO  N (0 ,  I) 

If  X  is  a  multivariate  normal  random  variable  with 
mean  vector,  y_,  and  covariance  matrix,  Z_  then 

Z  =  P  (X  -  y) 

is  multivariate  normally  distributed  with  mean  vector,  0_, 
and  covariance  matrix,  PEP '  [Ref.  8].   In  addition,  since  E_ 
is  the  symmetric  matrix  of  a  positive  definite  quadratic 
form,  there  exists  an  orthogonal  matrix,  B,  such  that 

BEB'  =  D 

where  D  is  a  diagonal  matrix  with  all  diagonal  elements 
positive  [Ref.  9].   The  matrix,  B,  can  be  constructed  using 
the  characteristic  vectors  of  Z_  as  columns  of  B.   Further, 
if  the  matrix  C  is  defined  as  the  diagonal  matrix  which  has 
diagonal  elements  equal  to  the  inverse  of  the  square-root 
of  the  corresponding  element  of  D,  then 

CBEB'C  =  CDC  =  I 
where  I_  is  the  identity  matrix.   Thus,  if  we  define 

P  =  CB 
then 
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Z  =  P(Y  -  y) 

is  multivariate  normal  with  zero  mean  vector  and  identity 
covariance  matrix,  and  the  elements  of  the  vector  Z^  are 
mutually  stochastically  independent  standard  normal  random 
variables . 

B.  KOLMOGOROV-SMIRNOV  TEST  OF  GOODNESS  OF  FIT 

Let  F(x)  be  defined  as  the  cumulative  distribution 
function  of  the  random  variable  X  which  is  N(0,1).   In 
addition,  define  S  (x)  to  be  the  sample  cumulative  distri- 
bution function  based  on  a  set  of  n  observations  of  a  random 
variable  assumed  to  be  N(Ofl).   Then 

S  (x)  =  k/n 

n        ' 

where  k  is  the  number  of  observations  in  the  sample  which 
are  less  than  or  equal  to  x.   Then,  the  Kolmogorov-Smirnov 
test  statistic  D  is  defined  as 

D  =  max|F(x)  -  S  (x) | . 

2\.  XI 

Observed  values  of  D  can  be  compared  with  its  tabled  distri- 
bution to  determine  the  acceptability  of  the  normal  hypo- 
thesis [Ref .  10] . 

C.  TEST  FOR  EQUALITY  OF  SEVERAL  COVARIANCE  MATRICES 
Suppose  Y, ,  Y?,...,Y   are  p-component  multivariate 

normal  random  variables  with  distribution  denoted  N  (y_.  ,%_•), 
i  =  l,2,...,q.   In  order  to  test  the  hypothesis 


39 


Z_.  =  Z_.  for  all  i,j  =  1,2,. ..,q, 

based  on  q  samples  of  size  N.  from  the  distribution  of  Y., 
i  =  1,2, ...,q,  let  the  following  quantities  be  defined: 


n.  =  N.  -  1,   i  =  1 , 2 , . . . , q ; 
li 


n  =   Z  n. ; 
i=l 

n . 

*i  =  ■/,  di,k  -  V(*i,k  -  ii>'    i  ■  1-2 g; 

k=l 


and 


q 

A  =   Z  A.  . 
i-l"1 


Then  the  test  statistic  is 


3  /v 

W  =  k[  Z  n. (log|zJ  -  log|£. |)] 
i=l  x  1 


where 


k  =  1  -r  1  -L--  I]  2P2  +  3P  -  l 

i=l  ni    n   6(P+1)  ((3-1) 

Z  =  A/n 


and 


Z.  =  Ai/ni    i  =  l,2,..,,q. 

Asymptotic  expansion  of  the  distribution  of  this'  test 
statistic  results  in  a  distribution  described  by  the 
following  probability  statement  for  W, 
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Pr(W<w)  =  Pr(Xf<w)  -  c[Pr(X^+4<w)  -  Pr(Xf<w)]  -  0(n  3) 


where 


q  1    1  2 

p(p+l)  [(p-D  (p+2)  (  I    -|--~)-6(q-l)  (1-k)^] 


i  =  ln.   n 
l 
c  = 


48k2 


f  =  k(q-l)   (p+l)p 

2 
and  Xf  denotes  a  chi-square  random  variable  with  f  degrees 

of  freedom  [Ref.  8]. 

This  test  was  applied  to  the  Air  Force  data  as  described 
in  III  C,  with  the  resulting  values  of  W  =  478.5,  f  =  252 
and  c  =  5.75.   The  test  statistic  is  extremely  significant 
at  the  .01  level  and  the  hypothesis  was  rejected. 

In  addition,  the  test  was  used  on  the  operational  data 
as  described  in  IV  C  2  where  now 


£  =  A/(n-q) 


and 


E±  =  Ai/(ni-l)    i  =  1,2,... ,q 


instead  of  as  defined  above.   The  results  for  this  test 
were  W  =  936.2,  f  =  728  and  c  =  14.9.   Once  again  the  test 
statistic  is  extremely  significant  at  the  .01  level  and  the 
hypothesis  was  rejected. 


41 


D.   TEST  FOR  INDEPENDENCE  AMONG  A  SET  OF  NORMAL  VARIATES 

If  the  p-component  vector  Y  has  a  multivariate  normal 
distribution  described  by  N  (y_,  £_)  ,  then  the  test  of  zero 
covariance , 

E[(Yi  -  M±.)  CYj  "  Uj)]  =  0, 

for  all  i  ^  j  and  i,j  =  l,2,...,p  is  equivalent  to  a  test 
of  the  stochastic  independence  of  the  components  of  Y. 
Suppose  a  sample,  Y,,Y_2,...,Y  /  °f  n  observations  from  the 
distribution  of  Y  is  obtained,  then  let  V  equal  the  deter- 
minant of  the  correlation  matrix,  R,  defined  by  equation  (2), 
The  asymptotic  expansion  of  the  distribution  of  V  results 
in  the  probability  statement  [Ref.  8]  i 


o  o  o  o 

Pr(-mlogVXv)  =  Pr(xf<v)  +  -y [Pr (xf+4^v)  -  Pr (xf^v) ]  +  0 (m~  ) 

m 

where 

2p  +  11 

m  =  n  -  -X— , 

6 

f  =  Jgp(p-l) 

and 

c  =  P(P-D  (2p2  -  2p  -  13) 

288 


This  test  was  applied  to  the  data  for  the  Pensacola  lab- 
oratory accumulated  in  the  Air  Force  experiment  as  des- 
cribed in  section  III  D.   Table  2  gives  the  numerical 


42 


results  for  the  ten  groups  of  samples  tested  and  the  test 
outcomes  for  a  level  of  significance  of  .10. 


TABLE  2 

RESULTS  OF  TEST  FOR  INDEPENDENCE 
f  =  21  m  =  5.833 

c  =  10.35        a  =  .10 


Sample  Group 

Test  Statistic 

Result 

1 

26.20 

Accept 

2 

39.10 

Reject 

3 

30.81 

Reject 

4 

31.57 

Reject 

5 

38.07 

Reject 

6 

36.36 

Reject 

7 

23.78 

Accept 

8 

41.60 

Reject 

9 

23.63 

Accept 

10 

35.17 

Reject 
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APPENDIX  B 

REGRESSION:   ESTIMATION,  TESTS  AND  PREDICTION 

A.   ESTIMATION  OF  REGRESSION  PARAMETERS 

Let  Y  be  a  seven-component  random  vector  with  distri- 
bution N (BX,  I) ,  where  B  is  an  unknown  7x2  matrix  of 
coefficients,  X  is  of  the  form  (l,x)  ',  and  £_  is  an  unknown 
covariance  matrix  which  is  constant  for  all  values  of  X. 
Then,  the  unbiased  estimate  of  B  [Ref.  8],  based  on  a 
sample  of  size  n  with  the  i —  from  N (BX . ,  E_)  is 

-1 
B  =  CA 

where 

n 
C  =   E   Y.X. 
"    i=i  -1"1 

and 

n 
A  =  I  x.x!  . 
i=l  -1"1 

This  estimate  is  normally  distributed  with  mean  matrix,  B, 
and  covariance  matrix  E_  ®  A   ,  where  the  symbol,  <8,  denotes 
the  Kronecker  product  [Ref.  11]  .   Further,  the  unbiased 
estimate  of  E_  [Ref.  8]  is 

1    n 
Z  =  -^      I       (Y.  -  BX.)(Y.  -  BX.  )  ' 
—    n-2  ._,   —l    — l   — l    — l 

and  (n-2)£  has  a  Wishart  distribution  with  parameters  E_ 
and  n-2. 
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B.  TEST  OF  THE  VALUE  OF  REGRESSION  COEFFICIENTS 
Suppose  the  7x2  matrix  B  is  partitioned  in  the  form 

B  =  (§!,  B2)  . 

Then  for  any  non-null  vector  of  constants/  C,  the  hypothesis 

C'B_2  =  C'B* 

can  be  tested  using  an  F  statistic  [Ref .  12] .   By  choosing 
C  to  be  the  vector  with  1  as  the  i —  component  and  all  other 
components  zero,  and  the  vector  B*  to  be  the  null  vector, 
the  hypotheses 

B2i  =  °     i  =  1/2,. .  .  ,7 

can  be  tested.   With  this  selection,  the  standard  t-test 
of  the  slope  of  a  regression  line  [Ref.  6  and  7]  can  be 
used.   In  this  way,  each  of  the  components  of  B^  can  be 
tested  individually,  each  with  an  assigned  level  of 
significance,  a. 

C.  CONSTRUCTION  OF  THE  REGION,  R  (x) 

a 

If  the  matrices  B  and  Z_  are  estimated  with  a  sample  of 
size  n  in  the  manner  described  above,  and  if  Y  is  a  new 
observation  from  N  (BX* ,  Z_)  ,  the  confidence  region  R  (x)  can 
be  constructed  so  that  Y  will  be  in  R  (x)  with  probability 
at  least  1-a.   The  estimate  of  the  mean  of  Y  is 


Y  =  BX* 
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and  has  covariance  matrix, 


where 


and 


-_1  -1 
S  =  T(A  8  £   )   T* 


n 
A  =   I   X.X! 
"    i=l  -1"1 


T  =  X*  8  I  . 


Thus, 


(Y  -  BX*) '  S  1(Y  -  BX*) 


2 

has  Hotelling's  T   distribution  [Refs.  8  and  5].   Hence,  the 

set  of  vectors  m  satisfying 

(Y  -  m)  •  S_1  (Y  -  m)  <  T2  (a) 


comprise  a  100  (l-a)%  confidence  region  for  BX*  [Refs.  8  and 
5].  Since  the  region,  R  (x) ,  is  to  place  bounds  on  Y,  which 
has  covariance  matrix,  £_,  it  can  be  defined  by 

(Y  -  m)  '  (S  +  l)        (Y  -  m)  £  T2  (a)  . 


The  set  of  all  vectors,  m,  satisfying  this  constraint  form 
a  confidence  region,  R  (x) ,  which  will  contain  Y  with 
probability  at  least  1  -  a  if  Y  came  from  N(BX*,  Z). 
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CCPPUTER    CUTPUT 


A.   DATA  FRCP 

AIRFCRCE 

EXPERIMENT 

GRCIP  A   1 

AL 

FE 

CR 

AG 

CU 

HG 

NI 

3 

76 

12 

4 

62 

■a 

5 

1 

12 

8 

3 

58 

3 

1 

C 

79 

6 

2 

63 

■a 

1 

1 

79 

10 

4 

61 

3 

5 

1 

73 

4 

2 

61 

2 

1 

0 

75 

9 

4 

61 

■a 

1 

3 

73 

11 

4 

59 

3 

2 

1 

79 

3 

2 

66 

3 

1 

0 

74 

2 

3 

62 

•a 

mi 

0 

GRCUP  k       2 

AL 

FE 

CR 

AG 

CU 

NG 

M 

C 

127 

7 

1 

12 

10 

2 

C 

142 

1C 

2 

13 

11 

3 

1 

144 

15 

3 

14 

11 

3 

0 

147 

9 

2 

14 

11 

3 

c 

142 

9 

3 

14 

11 

3 

c 

154 

14 

-a 

15 

12 

4 

c 

145 

17 

4 

15 

11 

4 

c 

153 

14 

3 

15 

13 

4 

2 

77 

11 

5 

6 

14 

1 

C 

84 

8 

4 

8 

16 

2 

GRCUP  *   3 

AL 

FE 

CR 

AG 

CU 

MG 

NI 

C 

151 

14 

3 

16 

12 

3 

1 

151 

13 

2 

14 

11 

4 

1 

84 

9 

4 

8 

16 

2 

c 

84 

12 

4 

6 

15 

1 

c 

ec 

10 

4 

c 

14 

2 

c 

82 

7 

4 

7 

15 

1 

c 

fl  ■* 

11 

5 

8 

16 

1 

c 

88 

6 

4 

7 

17 

0 

1 

77 

7 

4 

6 

14 

1 

1 

82 

8 

3 

8 

15 

1 

GRCLP  #   4 

• 

AL 

FE 

CR 

AG 

CU 

PC 

M 

26 

27 

9 

4 

35 

1C 

13 

23 

27 

8 

4 

26 

11 

13 

24 

27 

15 

4 

37 

1C 

14 

20 

28 

11 

4 

34 

11 

13 

24 

2e 

15 

5 

34 

1C 

14 

23 

29 

12 

5 

34 

11 

14 

22 

28 

10 

5 

24 

1C 

13 

23 

28 

14 

5 

35 

10 

14 

GRCLP  U       5 

AL 

FE 

CR 

AG 

CU 

*G 

M 

2 

137 

9 

4 

ICO 

6 

1 

5 

132 

17 

6 

107 

c 

1 

6 

127 

17 

6 

103 

5 

1 

3 

133 

2 

4 

116 

1 

A 

141 

4 

3 

111 

6 

0 

5 

-* 

138 

6 

4 

117 

6 

1 

2 

139 

12 

5 

1C3 

c 

1 

4 

145 

18 

6 

114 

6 

2 

3 

137 

11 

5 

1C9 

5 

1 

3 

132 

1C 

5 

1C4 

6 

1 
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GRCUP  * 

t 

AL 

FF 

CP 

AG 

CU 

*G 

MI 

r 

112 

11 

2 

13 

19 

4 

c 

105 

c 

2 

11 

19 

2 

1 

107 

9 

3 

13 

19 

2 

2 

gq 

10 

4 

13 

17 

3 

C 

IP  9 

8 

2 

12 

18 

2 

C 

104 

11 

3 

15 

lc 

3 

<», 

110 

1C 

3 

1? 

20 

4 

r 

115 

11 

3 

15 

21 

3 

1 

110 

13 

3 

1  i. 

20 

3 

1 

125 

11 

3 

16 

22 

3 

GRCUF  m 

- 

AL 

FE 

CR 

AG 

CU 

VG 

NI 

3 

95 

4 

4 

83 

27 

1 

t. 

95 

8 

A 

7Q 

25 

1 

1 

87 

10 

2 

72 

25 

2 

O 

.> 
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11 

4 

R6 

27 

1 

i 

1  Q4 

10 

4 

9C 

28 

1 

10  "' 

10 

5 

7Q 

27 

2 

■a 
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11 

5 

33 

2? 

5 

£ 
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11 

4 

37 

28 

1 

i 

104 

8 

4 

81 

27 
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GROUP  * 

8 

AL 

FF 

CR 

AG 

CU 

*G 

NI 

2 
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6 

3 

5 

11 

1 

c 
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8 

5 

5 

11 

3 

i* 

127 

4 

3 

8 

11 

3 

r 

133 

11 

4 

12 

4 

1 

143 

4 

2 

6 

12 

2 

1 

138 

1 

2 

6 

13 

1 

C 

127 

2 

3 

6 

11 

2 

r 

130 

2 

2 

8 

12 

2 

r 

137 

7 

4 

7 

12 

3 

r 

123 

Q 

4 

3 

11 

4 

GRCUP  * 

o 

AL 

FE 

CR 

AG 

CU 

MG 

NI 

2 

17 

12 

3 

9 

12 

1 

1 

la 

11 

2 

10 

12 

0 

■a 

17 

11 

2 

10 

13 

1 

c 

17 

3 

1 

10 

12 

0 

1 

17 

13 

2 

10 

12 

0 

1 

19 

1C 

2 

11 

13 

1 

i 

13 

6 

1 

11 

12 

1 

l 

18 

9 

2 

11 

12 

1 

CROUP  * 

10 

AL 

FE 

CP 

AG 

CU 

*G 

NI 

f 

137 

8 

3 

7 

10 

4 

1 

183 

10 

3 

8 

10 

7 

* 

197 

12 

3 

o 

9 

4 

r 

190 

8 

3 

7 

1A 

4 
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0 

n 

r, 

o 
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188 

8 

2 

6 

10 

4 

1 

212 

7 

2 

10 

1€ 

3 

f 

206 
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1 

10 

11 

4 

c 

217 

8 

2 

10 

11 

5 
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B.   SAMPLE  MEANS  ANE  STANCARC  CEVIATIGNS 


grcup  a 

MEAN 
S.D. 


I: 

AL 
l.ll 
1.  167 


EE 
75.56 

2.83* 


CR 
7.22 
3.632 


AG 
3.11 

C.928 


cu 

61.44 
2.299 


MG 
2.89 
0.333 


M 
1.89 
1.833 


GR01P  #       2: 
AL 
MEAN     0.20 
S.C.    0.422 


EF 
145.70 

7.8-8 


CR 
12.20 
3.225 


AG 
2.60 

0.843 


CU 
14.20 
1.135 


MG 

11.30 
0.823 


M 

3.30 
0.675 


GROLP  H       3: 

AL  Et  CR 

MEAN     C.50  82.10  8.90 

S.D.    C.7C7  3.384  2.C25 


AG 
4.10 

C.568 


CI 
6.90 
I  .101 


15.20 
1  .033 


M 
1.20 
0.632 


GRCLP  H       4: 

AL  FE 

MEAN    23.12  27.75 

S.D.    1.727  C.7C7 


CR 
11.75 
2.712 


AG 
4.50 
0.535 


CU 
34.88 
1.126 


MG 
1C.37 
0.518 


M 

13.50 
C.535 


GRCLP  H 

MEAN 
S.D. 


AL 
3.50 
1.269 


FE       CR       AC       CL       MG       M 
136.20    1C.6C     4.6G   108. 4C     5.50     l.OC 
5.119    5.582    1.033    5.969    0.527    C.471 


GRCLP  A   6: 
AL 
MEAN     0.5C 
S.D.    C.70  7 


FE 
1C9.70 
7.530 


CR 
9.9C 
2.183 


AG 
2.8C 
0.632 


CU 
13.  7C 
1.567 


MG 
19.40 
1.43C 


M 
2.90 
C.738 


GRCLP  U       7: 
AL 
MEAN     3.11 
S.C.    1.537 


FE 
ICO. 22 
6.039 


CR 
9.22 
2.279 


AG 
4.CC 
C  .866 


CU 
82.78 
5.653 


MG 
26.89 
1  .168 


NI 
1.67 
1.323 


GRCUP  * 

MEAN 
S.C. 


8: 
AL 
C.4C 
0.699 


FE 
131.40 
5.702 


CR 
5.4C 
3. 340 


AG 

3.20 
1  .C33 


CU 
6.6C 
1.174 


MG 
11.60 
C.699 


NI 
2.50 

1.C8C 


GRCUP  *   9: 

AL  FE 

MEAN     1.25  17.62 

S.C.    0.886  C.744 


CR 
10. CO 
2.268 


AC 
1.67 
C.641 


CU 
10.25 
0.7C7 


MG 
12.25 
0.463 


NI 

C.63 
.518 


GRCUP  k    10: 
AL 
MEAN     C.33 

S.C.    0.5  00 


FE 
198.56 
11.294 


CR 

8. CO 
2. 062 


AC 
2.33 

C.7C7 


CL 

8.78 
1.2C2 


MG 
10.11 

0.601 


M 
4.33 
1.118 
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C.   SAMPLE  CCRPELATICN  MATRICES 


GRCLP 

*       l: 

AL 

FE 

CR 

AG 

CL 

MG 

NI 

AL 

l.COO 

-C.210 

0.642 

0.449 

-0.254 

0.C36 

C.591 

FE 

-C21C 

l.COC 

-C.C50 

-C.169 

0.746 

0.338 

C.350 

CR 

0.642 

-C.050 

1.000 

0.808 

-0.5C7 

0.333 

C.736 

AG 

C.449 

-C.  169 

0.808 

l.CCC 

-0.495 

0.449 

C  .596 

CU 

-0.254 

C.746 

-0.5C7 

-0.495 

l.OCC 

0.073 

-C.C76 

MG 

C.C36 

0.338 

0.333 

0.449 

0.073 

1.000 

C.182 

NI 

0.591 

C.35C 

C.736 

0.596 

-0.076 

0.182 

1.000 

GRCLP 

H       2: 

AL 

FE 

CR 

AG 

CU 

MG 

NI 

AL 

1  .000 

C.121 

0.294 

-C.C62 

-0.093 

-C.192 

C.156 

FE 

0.121 

1.000 

0.631 

0.568 

0.831 

0.8C7 

C.816 

CR 

C294 

0.631 

l.OCC 

C.809 

0.747 

0.519 

0.735 

AG 

-0.062 

0.568 

C.8C9 

1.000 

0.789 

C.512 

C.625 

CL 

-C.093 

C.831 

0.747 

0.789 

1.000 

0.761 

C.636 

MG 

-C.192 

C.607 

C.519 

C.512 

0.761 

1.000 

C.620 

M 

0.156 

0.R16 

0.735 

C.625 

0.638 

C.620 

l.OCC 

GRCLP  H 


AL 


FE 


CR 


AG 


CL 


MG 


M 


AL 

l.CCC 

-C.627 

0.116 

0.138 

-0.C71 

-0.456 

0.000 

FE 

-0.627 

1.000 

-C.225 

-C.237 

C.510 

0.916 

-0.218 

CR 

C.116 

-0.225 

l.OCO 

C.493 

-0.254 

-C.308 

C.278 

AG 

C.138 

-C.237 

0.493 

l.OCO 

-0.160 

-0.C38 

-C.062 

CU 

-C.071 

C.51C 

-0.254 

-C.16C 

l.CCC 

C.704 

C.C32 

MG 

-0.456 

C.916 

-0.308 

-0.038 

0.7C4 

1.000 

-C.238 

NI 

C.COC 

-0.218 

0.278 

-0.062 

0.0  32 

-0.238 

1.000 

ROUP 

*       A: 

AL 

FE 

CR 

AG 

CU 

MG 

NI 

AL 

1.000 

-C.439 

0.C69 

-C.C77 

0.377 

-C.54C 

C.232 

FE 

-C.439 

l.OOC 

0.261 

0.756 

-0.763 

0.293 

C.378 

CR 

C.C69 

C.261 

l.CCO 

C.394 

0.082 

-0.433 

C.887 

AG 

-0.077 

C.756 

0.394 

l.CCC 

-0.593 

-0.258 

C.50C 

CU 

C.377 

-C.763 

0.082 

-0.593 

1.000 

-0.153 

0.119 

MG 

-0.54C 

0.293 

-0.433 

-C.258 

-0.153 

l.COC 

-C.258 

M 

C.232 

0.378 

0.887 

C.50C 

0.119 

-0.258 

l.COO 

GRCLP  U      5: 


AL 


FE 


CR 


AG 


CL 


MG 


M 


AL 

l.CCC 

-C.444 

0.518 

0.5C9 

0.044 

-0.249 

c.ooo 

FE 

-C.444 

1.000 

-0.  110 

-0.286 

0.369 

C.536 

0.184 

CR 

0.518 

-cue 

1.000 

C.929 

-0.385 

-0.227 

C591 

AG 

C.5C9 

-0.286 

0.929 

l.OCO 

-0.256 

-0.408 

0.685 

CU 

C.044 

C.369 

-0.38  5 

-C.256 

l.OCO 

0.142 

C.119 

MG 

-C.249 

0.536 

-0.227 

-0.408 

0.142 

l.OOC 

0.000 

NI 

C.OCC 

C  .  1  84 

0.591 

0.685 

0.119 

0.000 

1.000 
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GROUP  «   6: 


AL 


Ff 


CP 


AG 


CU 


MG 


MI 


AL 

l.ooo 

-C.19P 

0.252 

j.745 

3  •  0  5  0 

-C.220 

-0.106 

FE 

-J. 198 

l.roc 

0.322 

-0,200 

^.575 

0.890 

0.214 

CF 

.1.252 

C.  322 

l.GOO 

0.467 

1,737 

0.37C 

".614 

AG 

3.7^5 

-C.2CC 

0.467 

1.000 

0.^-93 

-0.02^ 

0.191 

CU 

0.0  50 

r.F75 

0.737 

0.493 

1.000 

0.7C4 

0.548 

MG 

-0.220 

Q.89C 

0.370 

-C.024 

3.7C4 

L.000 

1.253 

NI 

-0.106 

c.214 

0.614 

0.191 

T.^48 

0.253 

1.000 

group  n 


AL 


FE 


CR 


AG 


CU 


mg 


NI 


AL 

l.OfO 

C  .  4  6  9 

0.206 

0.376 

^.593 

C  .426 

-1.164 

FE 

3.469 

1.  )0C 

0.2C5 

0.7  +  1 

).?64 

0.90° 

"».^RQ 

CP 

.3.206 

t.2^5 

1.000 

0.053 

-1.122 

1.246 

3.359 

AG 

1.376 

0.741 

0.C63 

1.000 

^.^60 

0.618 

3.32^ 

CU 

:.593 

^.76^ 

-0. 122 

0.460 

i.  ono 

0.773 

-:.228 

MG 

J. 426 

C  .  9<  9 

0.2^6 

o  .618 

1.773 

1.000 

1.216 

NI 

-0.164 

CO  89 

n.359 

1.327 

-3.228 

0.216 

l.oon 

GR^uP  tf       ^: 


AL 


FE 


CR 


AG 


CU 


MG 


NI 


AL 
FE 
CP 
AG 
CU 
MG 
NI 


1.000 
1.457 
-3.219 
-3.431 
-0.5C6 
3.136 
-3.7?b 


C.^c7 
l.C  *  C 
■C.196 
•C.5^-3 
•0«G*G 
C.770 
■0.343 


-0.219   - 


196 
GOO 
812 
045 
3C4 
0.770 


0. 

1. 

0. 
y. 
0. 


431 
543 
312 
000 
0.110 
0.492 
n.697 


1 


-0.596 

-1.040 

1.045 

-0.110 

l.OOC 
5.054 
3.5  26 


0.136 
0.770 
■0.304 
•0.492 
0.O54 
1.000 
■3.294 


.736 
.343 


u.77C 

3.697 
0.526 
0.294 

1.0  JO 


GROUP    H 


9: 
AL 


Ft 


CP 


AG 


CU 


MG 


MI 


AL 
PE 

CP 
*G 
CU 
MG 
NI 


UGOO 

-3.271 

1.^26 

3.566 

-3.342 
3. =22 
3.5*t5 


-(  .  271 
1 .  •'  C  0 
■0.339 
■C.112 
0.747 
0.311 
C325 


1.426 

0.566 

-0, 342 

0.522 

1.545 

0.3  39 

-C.112 

).747 

0.311 

0.325 

l.  000 

0.736 

-1.624 

0.136 

-0.243 

0.  7  36 

1.000 

-0.552 

0.12C 

0.269 

3.624 

-0.5  52 

1.000 

3.218 

0,293 

3.  136 

0.120 

0.218 

1.000 

0.447 

0.243 

0.269 

3.293 

0.447 

l.OOC 

GROUP    ft    10: 


AL 


r-F 


CR 


AG 


CU 


MG 


NI 


AL 

1.000 

C.03G 

0.606 

0.354 

0.139 

-0.555 

^.224 

FE 

1.C3C 

1.    00 

-0.376 

-C.621 

).830 

1.506 

-3.264 

CR 

3.606 

-C.376 

l.COO 

0.772 

-0.303 

-0.706 

0.38C 

AG 

3.354 

-C.621 

0.772 

l.OOC 

-3.734 

-0.686 

1.316 

CU 

3.139 

C.83C 

-0.303 

-0.784 

1.0  00 

0.38  5 

-1.217 

MG 

-1.555 

r.5C6 

-0.7C6 

-0.6  86 

0.38  = 

1.000 

1.124 

NI 

0.224 

-C.264 

0.380 

0.316 

-0.217 

0.124 

1.000 
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D.   RESULTS  CF  REGRESSION  CF  VARIANCE  CN  PEAN  SQUARED 


ALLVINUN! 

INTERCEPT: 

SLOPE: 

CORRELATION: 


0.8768E  OC 
0.4C16E-C2 
C.7C67E  CC 


RESICUAL  VARIANCE 

VARIANCE  OF  INTERCEPT 

VARIANCE  OF  SLCPE 


C.AU8E  CI 
0.5782E-C1 
C.1824E-C5 


IRON 

INTERCEPT 
SLOPE: 
CORPELATICN 


-0.4542E  CI 

0  .3C8CE-C2 
0.9296E  OC 


RESIDUAL  VARIANCE 

VARIANCE  OF  INTERCEPT 

VARIANCE  CF  SLCPE 


0.1829E  C4 
C.5537E  C2 
C.83C5E-C7 


CFRCf  ILM 

INTERCEPT 
SLCPE 

COPRELATICN 


C.7276E  01 
0  .2642E-C1 
0.1193E  OC 


RESICUAL  VARIANCE 

VARIANCE  OF  INTERCEPT 

VARIANCE  CF  SLCPE 


C.5985E  C3 
C.5717E  C2 

C.8870E-03 


SILVER 

INTERCEPT 

SLOPE 

COPRELATICN 


0.5494E  OC 
0.7357E-C2 
0.1654E  CC 


RESICUAL  VARIANCE 

VARIANCE  OF  INTERCEPT 

VARIANCE  OF  SLCPE 


0.7648E  00 
C.A397E-01 
C.58C0E-04 


CCPPEP 

INTERCEPT 

SLOPE 

COPRELATICN 


0.3822E  CC 
0.3238E-C2 
0.9494E  OC 


RESIDUAL  VARIANCE 

VARIANCE  OF  INTERCEPT 

VARIANCE  OF  SLCPE 


C.1631E  03 
C.2883E  CI 
0.1054E-C6 


PAGNES. 

INTERCEPT: 

SLOPE: 

COPRELATICN: 


0.2444E  CC 
0.2224E-C2 
0.7544E  CC 


RESICUAL  VARIANCE 

VARIANCE  CF  INTERCEPT 

VARIANCE  CF  SLCPE 


C.15C3E  CI 
0.373AE-01 
C.2511E-C6 


NICKEL 

INTERCEPT 

SLOPE 

CORRELATICN 


0.1C66E  CI 
-0  .4177E-C2 
-0.237CE  OC 


RESICUAL  VARIANCE 

VARIANCE  OF  INTERCEPT 

VARIANCE  OF  SLCPE 


C.8284E  CI 
C.1239E  00 
C.3126E-C4 


RESULTS  OF 
REGRESSION 


A  T  TEST  CF  THE  HYPOTHESIS 
LINE  IS  ECUAL  TC  ZERO 


THAT  THE  SLCPE  OF  THE 


ALUMNUN 

IRCN 

CHRCMIUN 

SILVER 

COPPER 

MAGNES. 

NICKEL 


0.29735A2E  01 
0.1C68822E  C2 
0.8869581E  00 
C.9660166E  00 
0.9976842E  01 
0.4439077E  01 
0.7471411E  00 


REJECT 
REJECT 
ACCEPT 
ACCEPT 
REJECT 
REJECT 
ACCEPT 
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E.   OPERATIONAL  CATA,  ESTIMATED  PARAMETERS  AND  TESTS  OF 
REGRESSION  SLOPES  FOR  ENGINES  CF  THE  R182C-82  NCCEL 

RESULTS  CF  LINEAR  REGRESSION  USING  THE  NCOEL, 

Y  =  BX  ♦  E 
WHERE: 

Y  IS  A  VECTCR  CF  WEAR  METAL  CCNT AM  I  NAT  I CN  IN  CIL 

SANPLES 

Yl  =  PPM  CF  ALUMINUM 

Y2  =  PPM  CF  IPCN 

Y2  =  PPM  CF  CFRCMIUM 

>4  =  PPM  CF  SILVER 

Y5  =  PPP  CF  CCFPER 

Vt  =  PPM  CF  MAGNESIUM 

Y7  =  PPM  CF  MCKEL 

B  IS  A  7X2  MATRIX  CF  COEFFICIENTS 

X  IS  A  VECTCR  CF  VARIABLES 

XI  =  1 

X2  =  HCURS  SINCE  OIL  CHANGE 

ANC  t  IS  U  CCMPCNENT  MULTIVARIATE  NORMAL  RANDOM 
ERRCR  VECTOR  UITH  MEAN  VECTCR  2ERC  AfsD  UNKNOWN 
COVARIAfsCE  NATRIX  SIGMA 
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DATA 


ENGINE    NU^B1 

5  R :     5 

1672V 

• 

'yi 

Y2 

Y3 

Y4 

Y5 

Y6 

Y7 

XI 

6.- 

25.0 

&.0 

2.0 

12.' 

3.C 

1,0 

11.0 

£.r 

28.0 

8.0 

2.0 

11." 

1.0 

2.C 

57.  1 

10. 

41.0 

3.  1 

3.-» 

2  3.r 

2.C 

2.C 

256.  : 

9. 

25.0 

1.0 

2.3 

12.0 

}.: 

1.0 

21.0 

s.r 

3 1 . 1 

6  .  0 

2.0 

10.-, 

2.o 

l.C 

OA.  ^ 

K. 

30.0 

c.r 

1.0 

16.'^ 

2.0 

%r 

99.  ~ 

6.  . 

26.0 

3.  ' 

2.") 

14. ~ 

1.0 

w         0 

64.0 

11.'- 

?  ).n 

1  .0 

1.0 

16. <^ 

1.0 

l.f 

133.  n 

7.. 

27.? 

3.'~> 

12.- 

l.C 

2.^ 

108.0 

11. 

3  8,0 

7.0 

1. 1 

19.  .0 

2." 

?.C 

171." 

NUVBCR  OF  p^TA  PCIMS  =  1 


RESULTS 


PI 


ESTIMATE  C  e  R 
P2 


7.14706 

'.  .:  15  62 

23.26596 

0.06875 

A  ,  '  -7  C,  P  c 

-0.00  24  1 

1.  5C8(  A 

".COO  Q3 

9.7187C 

?. . :  4  e  l  c 

1.5472? 

0.00053 

1.93402 

!.C0n66 

ESTIMATE  CF  CCVARIANCC  VATRlx  SIGMA 


2.5<2 

1.314 

-1.C98 

-0.237 

.)  .  4  r'  9 

-0,  ->75 

-w.9^7 

1.  314 

4.626 

3.390 

-^.792 

"■  .  8  8  3 

0  .  ^  R  c 

-r.65^ 

1 .  C  9  8 

3.393 

8.665 

'..5  3  8 

-2.176 

6.658 

-0.615 

C.237 

J.7Q2 

0.538 

^.705 

1.482 

".172 

-^.004 

CATC 

0.883 

-2.176 

0.432 

4,619 

0.4-71 

e . q :  8 

C.^-75 

^.955 

(..65  8 

0.  172 

0.471 

^.s^s 

-C.252 

C.93"7 

-C.650 

-".615 

-  U094 

'.poo 

-0.252 

1.74  7 

T  TEST  rr-     THE  HYPOTHESIS  THAT  4 
COMPONENT  CF  THE  62  VCCTCP  EQUALS 
ZE°C  AT  A  .10  SIGNIFICANCE  LEVEL 


TEST  STATISTIC 

2.15173 
7.C1834 
-0. 17979 
T. 22784 
4.c1423 
C. 15746 
r.iir25 


RESULT 

REJECT 
RCJECT 
ACCEPT 
ACCEPT 
REJECT 
ACCEPT 
ACCEPT 
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ENGINE  NIJWP6P:  515355 


C&TA 


Yl 


Y2 


Y3 


Y^ 


Y5 


Y6 


Y7 


MJ^PrF    Df    CAT4    POINTS    = 


XI 


K  . 

3  9 .  G 

5.0 

0,0 

17.0 

2.0 

4.C 

126." 

9.r 

2P.0 

2.0 

l.o 

1 1 .  0 

2.C 

3.n 

117,0 

12.* 

3 1 .  G 

3.0 

Q.T 

14.  ~ 

2.C 

2.C 

84.  0 

6.0 

23.0 

f.O 

2.0 

10.0 

o . ; 

'^.0 

66.0 

5«i 

19.0 

<3.n 

c.  ■) 

7  . ' , 

0.0 

l.C 

24.0 

4. 

2  J.r 

3.0 

1.  ^ 

7.0 

1.0 

2.C 

10,0 

P. 

33.0 

l.o 

0.  3 

15. n 

2 .  : 

1.0 

7Q.  r 

-  •  . 

26.  J 

6 .  n 

2.  3 

1 2 .  o 

1.0 

i.e. 

46.0 

lt.f 

2  8.0 

g.n 

3.  1 

11.0 

2.0 

4.0 

l<5.r 

P f 5ULTS: 

FS 

flMATE    CF    E 

PI 

H2 

K. 19364 

G. 0  38 9  3 

lc. H3*  9C 

0.11819 

6.56766 

-0.03657 

1.73114 

-0.01152 

7.73^36 

0.06015 

0.£P'42 

0. 01023 

1.43462 

0.0C851 

re 

TIMATE    CF    CCVARIANCc    MATRIX    SIGMA 

5.787 

6..  406         2.578 

0.336          2.786 

1.337 

2.013 

6.4f  6 

20.919         4.305 

0.46  6       10.39  7 

2.700 

3.^44 

2.5  78 

4.3'^5         9.^61 

-0.004          1.62C 

0,526 

2.66  5 

C.336 

0.466       -0.004 

1.160         1.115 

3.095 

C.35C 

2.7P6 

1J.397          1.620 

0.115          5.864 

1.090 

^.774 

1.337 

2.70",          0.526 

O.095          1.090 

0.646 

0.816 

2.r13 

3.444          2.665 

0.350         0.774 

0.816 

2.125 

T  TEST  CF  THE  hYPOTHESIS  THAT  A 
COMPONENT  CF  THE  62  VECTCP  EQUALS 
ZEPC  AT  A  .K  SIGNIFICANCE  LEVEL 


TEST    STATIST! 

1.92706 

3. r7322 
-1.  <^  2^41 
-1.27241 

2.C541A 
1.51371 
0.72-d72 


RESULT 


RE- J 
REJ 
ACC 
4CC 
REJ 
ACC 
ACC 


ECT 
ECT 
EPT 
EPT 
ECT 
EPT 
EPT 
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OATA 


Yl 


Y2 


ENGINE  NUMBER:  51667a 


Y3 


Y4 


Y5 


Y6 


MIN'BCF    OF    TATA    FCINTS    = 


Y7 


XI 


9.0 

25.0 

5.0 

1.0 

13.0 

3.0 

6.0 

60. o 

P.-' 

25.0 

7.0 

1.0 

1  ?.o 

l.C 

5.^ 

8<?.0 

F  .  ■"' 

25.0 

3.0 

1.0 

11.0 

1.0 

1.0 

79..; 

7.. 

2%-) 

2.0 

0,  3 

14," 

•2.~ 

5.0 

58.0 

1 C  .  1 

3  5.0 

5.0 

0.0 

14.0 

3.r 

3.0 

29.0 

P." 

31.C 

1.0 

0.0 

12.  r 

3.0 

4.0 

31.^ 

6.  0 

3  ) .  0 

8.C 

1.0 

12.0 

4.0 

^.<", 

117.0 

4.0 

25.  C 

5.C 

2.0 

m.n 

2.0 

3.0 

51. C 

2  3.0 

9.0 

1.0 

e.: 

1.0 

O.C 

18.  Q 

SESUt  T5: 


01 


ESTIMATE    OF    R 


6.  732  '7 
26.2383C 
3.a'  117 
0.4037^ 
1  1.05223 
1.85873 
2.' 8  8  7r 


n 


0 
r. 


P2 

j  1 2 : 5 

0  0  c  3  7 
017Q1 
-_:063  3 
C1226 
0  0  615 

,01429 


ESTIMATE    CF    CCVARIANCE    MATRIX    SIG^A 


i 

A 
-1 

-G 

2 


8  64 
322 

961 
956 
Cc3'- 
1.468 


6.  322 

27,427 

-2.397 

-1.02  3 

4.363 

3.946 

].959 


-1.535 
-2. 3^7 
7.<U5 
3.869 
-3.111 
-0.413 
-2.296 


-i 


961 
023 

860 

0.46  2 
010 
267 
438 


2.9  56 
4.368 
-3.111 
- 1 .  0  1  C 
4  .  Q  >-■  8 
1.119 
2,893 


_~« 


3.94A 
.413 
:,267 
1.119 
1.321 
0,80  3 


1.468 

r.95S 

-2.296 

-0.4  38 
2.893 

0.803 
4 .  Q  5  r 


T  TEST  OF  THE  HYPOTHESIS  THAT  A 
COMPONENT  CF  THE  62  VECTCP  ECUALS 
ZfPC  AT  A  .10  SIGNIFICANCE  LEV^L 


TEST  STATISTIC 

r. 55162 

C  .  r  t  ?  2  0 

0.57286 
f . 83792 
0.54810 
f% .  4  3 1 3  1 
0.63835 


RESULT 

ACCEPT 
ACCENT 
ACCEPT 
ACCEPT 
ACCEPT 
ACCEPT 
ACCEPT 
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ENGINE    NUMRcp:     515058 


DATA: 


Yl 


Y2 


Y3 


Y4 


Y5 


Y6 


Y7 


XI 


10.:.' 

23.0 

7., 

25.0 

4.r 

20.  ^> 

17. 

3  3.0 

6.0 

24.0 

7*1 

20.0 

6.r 

15.0 

26.0 

5.0 

13*0 

7.0 

1.  5 

8.0 

0.0 

2.C 

0  .  ) 

10.0 

2.0 

6.0 

0.0 

8.0 

1.0 

C.Q 

n.o 

7.'? 

2.0 

5.0 

1.0 

9.0 

1.0 

7.0 

1.? 

9.0 

0  .  G 

3.0 

1. 1 

6.0 

l.C 

8.0 

").  1 

7.0 

O.O 

3.0 

1.0 

4.0 

0.0 

2.0 

2  6.0 

3.^ 

13c.O 

2.0 

R3.0 

1..° 

73.^ 

l.C 

107.^ 

0  .  f. 

<s«.o 

3.C 

4^.^ 

l.C 

56.0 

n.n 

8  8.  0 

MJ^tR    PF    CATA    PCINTS    =       9 


RESULTS 


01 


£STIM£TC    CF    P 


P.c3^6' 

21  .2^82c 

7.14537 

O.opr  n 

6.1622c 

-£.290*7 

1.2137' 


P2 

-O.C  1.984 

).C32  2~, 
-0.  034-48 
-0.0056  5 
0  .  0  1  Q  5  5 
0.01*22 
0. 00307 


ESTIMATE    CF    CCVARIANCF    ^AT^IX    SIG^a 


17, 5f 9 
18. 182 
-7.36c 

-0.6CC 
C.686 
2,f73 

-C.173 


18.182 


3'"> 
-3 
-1 

1 


05c 
"02 
663 
58S 
725 
941 


-7.369 
-3.002 

6.738 

2.0  24 
-1.333 
-0.-605 


-0.6  09 
-1.66  3 
%059 
%275 
-0.115 
-0.163 
-0.152 


0.8  86 
2.588 
2.H24 
•■'•.115 
3.289 
0.237 
0.6C7 


2, 


1 

-1 


"73 
725 
333 

163 
^.237 

0.525 
0.^98 


-1..  • 
-0. 

'  . 

-0. 
c. 

1. 


173 
941 

60  5 
152 
607 
40  p 

448 


T  TEST  OF  THE  HYPOTHESIS  THAT  A 
CCVPONENT  CF  THE  B2  VFCTCR  ECUALS 
ZFRC  AT  4  .10  SIGNICICANCE  LEVEL 


TEST  STATISTIC 

-C.^5716 

r. 56631 

-1.2  8071 

-1.03937 
C.c8615 
1.89257 
r. 2^613 


RESULT 

ACCEPT 
ACCEPT 
ACCEPT 
ACCEPT 
ACCEPT 
ACCEPT 
ACCEPT 


57 


ENGINE:  NUMBER:  509090, 


TATA 


Yl 


Y2 


Y3 


2 1 .  C 

29.0 

3.0 

1ft.' 

2  2.0 

2.C 

IP." 

4-  1.  0 

o .  n 

17. 

38.0 

o.c 

9.  S 

2  9.0 

5.C 

e.o 

21.0 

7.0 

14.  ( 

37.0 

6.0 

e.c 

25. C 

0.0 

15.' 

4  4.0 

5.C 

12.^ 

33.0 

3,r 

!ft. 

33.0 

4.c 

Y4 


NUMPEF  LF  DATA  PCINTS  =  11 


RESULTS: 

EST  I 

MATE  CF  P 

PI 

P2 

9.01636 

O.C 3755 

18.23833 

n. 10673 

3.P^371 

-0.03490 

r. 73897 

-0.00351 

4.  1188* 

C. 041 5  7 

C. 37839 

&.*.">  60  5 

1.30356 

-0.00160 

Y5 


Y6 


Y7 


XI 


0.0 

3  .  0 

3.0 

l.C 

111.0 

1.0 

7.c 

2.0 

1 .  c 

90.0 

1.0 

12.0 

1.0 

1.0 

16  l.C 

j.0 

11.0 

l.C 

1.0 

183.0 

0.0 

7.0 

o.c 

l.C 

9C.C 

1.0 

7.* 

i.C 

1.0 

30.0 

0.0 

12.r; 

1.0 

:■ . n 

198. 0 

0.0 

6.n 

0.0 

l.C 

60.0 

1.0 

11.0 

2.0 

4.0 

156." 

0.0 

12.0 

l.C 

C  .r 

193.0 

0.  0 

13." 

2  .  0 

l.C 

188.0 

FSTIMATE  CF  CCVARI4NCE  MATRIX  SIGMA 


14«99<- 

r.  956 
-3.  165 

r.176 

-r,36e 
3.122 

C.78F 


1.956 
20.278 
-0.15  5 

0.297 

1.499 
-0. 153 

3.944 


-3.165 

-0.155 

6.752 

0.h27 

0.771 

-C.035 

.282 


0.176 
°.297 
".42  7 
0.195 
0.130 
'.133 
0 .281 


-0.366 
1.49  9 
0.771 
0.130 
0.939 

-0.114 
0.187 


3.122 
-r.153 
-3.035 

0  .  1 3  3 
-0.114 

0.929 

0  .  3  c  1 


v. 
0. 

■ 
X. 


789 

282 
281 
187 
351 

20  2 


T  TEST  CF  tHF  HYPOTHESIS  THAT  A 
COMPONENT  CF  THE  C2  VFCTCR  FCUALS 
ZFRC  AT  A  .10  SIGNIFICANCE  LEVEL 


TFST    ST\TISTIC 

1.80834 

4.^2040 

-r. 35791 

-1.' 3459 

8.00668 

1.17122 

-C. 27252 


RESULT 

ACCEPT 
PEJECT 
ACCEPT 
ACCEPT 
REJECT 
ACCEPT 
ACCEPT 


58 


DATA: 


Yl 


Y2 


ENGINE  NUMBER:  524739 


Y3 


Y4 


Y5 


Y6 


NUMBER  OF  C4TA  POINTS  =  14 


Y7 


XI 


16.0 

?2.0 

8.C 

2.0 

4.0 

0.0 

2.0 

2  8.0 

13.  ,3 

37.  3 

C.C 

o.^ 

13.0 

2.0 

3.*~ 

181.0 

16.  n 

35.0 

1.0 

2.0 

11.0 

5.n 

6.0 

32.0 

6.0 

29.0 

0.0 

3.0 

9.0 

2.0 

3.0 

158.  r 

5.0 

27.0 

1.0 

o.o 

7.0 

O.C 

l.C 

15.C 

16.0 

39.0 

1.0 

2.0 

ll.r 

5.0 

3.0 

6  2.  r 

16.  r 

^5.0 

2.0 

4.^ 

14.  C 

4.0 

1.0 

169.0 

13.0 

2  2.0 

2.^ 

2.0 

7.0 

2.0 

3.0 

6.0 

m.o 

31.  J 

C.C 

0.0 

6.n 

3.0 

3.0 

19. C 

17.  r 

3  9.0 

4.C 

3.0 

13.0 

3.0 

1.0 

169.0 

12. r 

37.0 

1.0 

0.0 

12.0 

3.0 

3.G 

102.0 

11." 

35. C 

2.0 

0.0 

11.0 

2.0 

2.0 

117.  n 

9.C 

34.0 

l.C 

0.0 

\c '.r. 

2.0 

3.r 

133.0 

13.0 

3  0.0 

13.0 

4.0 

14.0 

3  .  0 

7.C 

148.0 

RESULTS 


Bl 


ESTIMATE  EF  P 


11.92C93 
27.11922 

?.^5ee2 

0.895  75 
6.64235 
2.22275 
3.-^211 


P2 

0.C0456 
0.06821 
0.005*4 
0. 00  70  6 
0.03660 
0.00365 
-0.00119 


ESTIMATE  CF  CGVARIANCE  MATRIX  SIGMA 


15.669 
9.095 

6.728 
2.776 
3.56F 
3.183 
L.805 


9.095 
26.769 
-1.43^ 
0.710 
8.839 
6.396 
1.686 


L 

-1 

14 

r 
-0 


728 
434 
312 

104 
37  9 
975 


2.828 


2.776 
0.710 
3.  104 
2.383 
0.515 
0.7^7 
0.671 


3.568 
8.839 
0.379 
0.515 
6.309 
2.357 
1.551 


3.183 
6.3°6 
-0.O75 
0.747 
2.357 
2.390 
1.151 


0  .  « -•>  c 
1.686 
2.  820 
n.671 
1.551 
1.151 
3.237 


T  TEST  CF  THE  HYPGTHESIS  THAT  A 
COMPONENT  CF  THE  P2  VECTOR  ECUALS 
ZERC  AT  4  .10  SIGNIFICANCE  LEVEL 


TEST  STATISTIC 

0.27445 
3.1402  6 
C. 34273 
1.09012 
4. 1^994 
r. 56178 
-r. 15716 


RESULT 

ACCEPT 
REJECT 
ACCEPT 
ACCEPT 
REJECT 
ACCEPT 
ACCEPT 


59 


ENGINE    NIJMRER:     516309 


CATA 


Yl 


Y2 


Y3 


YA 


Y5 


Y6 


NUMBER  IF  CATA  PCINTS  - 


RESULTS: 


81 


ESTIMATE  CF  R 

82 


2.44113 

0,C75^6 

F. 17f 8  8 

".34263 

?.<  328c 

3.0C355 

1. ]  I960 

-  0  .  C 1 2  6  6 

4.7G6CC 

Cn7843 

1.148?^ 

0.03404 

(.74821 

0.00472 

Y^ 


XI 


7.< 

31.0 

CO 

0.0 

10.  C 

0.0 

3.0 

63.0 

7.0 

2  ).0 

8.0 

1.0 

8.0 

0  .  c 

2.0 

32.  0 

6.0 

19.0 

2.0 

0.0 

9.0 

0.0 

1.0 

26.0 

6.0 

21.0 

0.0 

3.0 

9.0 

0.0 

0.0 

37. ^ 

5.f 

31.0 

CO 

0.0 

7.0 

1.0 

1.0 

64.0 

5.C 

26.0 

6.0 

0.0 

9.0 

1.0 

2.0 

9  CO 

15.0 

5  3.0 

5.r' 

0.0 

17.0 

4.C 

0.0 

llc.^ 

3.'- 

23.0 

8.0 

2.0 

7.0 

0  .  0 

CO 

17.C 

4.C 

12.0 

CO 

1.0 

4.0 

0.0 

".0 

3  2.". 

ESTIMATE  CF  CCVARIANCC  ^aTRIX  SIGMA 


6.41? 
13. 140 

1.098 
-CC17 

5.261 

1.  152 
-C7*7 


13.14^ 

49.28* 

3.719 

13.528 

3.445 

-2.533 


1.C.Q8 

3.719 

14.206 

1.6*6 

2.3«4 
0.7OS 

0.121 


-{  ."17 
0.820 

1.646 

-  ,39  5 

-0.073 

C180 

-C2?8 


5.261 
13.528 

2.3°4 
-0.073 

6.123 

C9°8 
-0.482 


1.152 

3.445 

r.795 

CISC 
C  ^r  p 

0.493 

-0.638 


747 
533 
121 

2:8 

432 
63^ 
1.40C 


v    i 

—  O 


_  f> 


T    TEST    CF    THE    HYPOTHESIS    THAT    A 
C.rwPONFNT    OF     THE    62    VECTOR    ECUALS 
ZERC    AT     A    ,1C     SIGNIFICANCE    LEVFL 


TEST    STATISTIC 

2. £2637 
4.6^781 
CC3989 
-l.c2397 
3.r 2476 
A. 62768 
C  ?  8  0  8  6 


RESULT 

REJECT 
REJECT 
ACCEPT 
REJECT 
REJECT 
REJECT 
ACCEPT 


60 


PATA 


ENGINE  NUMBER:  516094 


Yl 

Y2 

Y3 

Y4 

Y5 

Y6 

Y7 

XI 

5.  0 

16.0 

3.C 

0.0 

10. C 

1.0 

1.0 

le.o 

10.  c 

34.0 

7.n 

1.0 

11.0 

2.0 

3.0 

94." 

7.: 

20.0 

6.0 

2.  1 

3.0 

1.0 

0.0 

62.  C 

6.<" 

23.0 

9,0 

2.0 

8.^ 

2.0 

^.0 

93.° 

7.{ 

21.0 

l.C 

0.0 

8.C 

1.0 

l.C 

58. n 

8.  . 

32.0 

9.C 

0.0 

10. c 

1.0 

3.0 

38.  * 

1 1 .  0 

3  ).0 

9.0 

3.0 

7." 

1.0 

2.0 

84.  rj 

4.< 

27.0 

7.0 

1.0 

16.0 

2.0 

6.C 

84.  0 

7. 

24.0 

CO 

0.0 

5.0 

1.0 

1.0 

5^.0 

NUN'PLR  OF  rATA  POINTS  = 


PESLLTS: 


pl 


ESTIMATE  CF  P 


4.62017 
12.3°3C2 
-1.43774 
-r. 2^837 
7.1661P 
0.49376 
-1. 03911 


R2 


.  C  3  6  8  3 
.20229 
00912 
0.01297 
0.02886 
0. Oil 90 
C04780 


j 


FSTIVATE  CF  CCVARIANCE  "ATRIX  SIC^A 


4.682 

3.21  J 

0.237 

-0.817 

-A. 107 

-0.551 

-2.*94 

3.21C 

6.565 

1  .  P  1  3 

-0.012 

-1.4C2 

-0.429 

-2.173 

0.237 

1.813 

6.177 

0.322 

2.375 

-r.078 

C.530 

C  .  8 1 7 

-0.012 

p.32  2 

0.737 

0.257 

n.176 

-0.013 

4.107 

-1.492 

2.375 

0.  257 

10.486 

0.80  3 

3.Q22 

r.551 

-0.^29 

-0.078 

0.176 

0.8^-3 

0.185 

0.45  2 

2.^04 

-2.173 

C.530 

-0.C13 

3.922 

0.452 

2.373 

T  TrST  CF  THE  HYPOTHESIS  TH»T  A 
COMPONENT  Cc  THE  R2  VECTCR  ECUALS 
ZFRC  AT  A  .10  SIGNIFICANCE  LEVEL 


TEST  STATISTIC 

1.203  53 
5.57493 
2.81619 
1.C665C 
C.t?926 
l.c5434 
2.19124 


RESULT 

ACCEPT 
REJECT 
REJECT 
ACCEPT 
ACCEPT 
REJECT 
REJECT 


61 


ENGINE    NIJMRcr:    516949 


CATA 


Yl 


Y2 


Y3 


Y4 


Y5 


Y6 


MJMPEF    OF    TATA    FCIMS    = 


i  o 


Y7 


XI 


5.i 

22.0 

1.0 

1.0 

io. r 

2.C 

5.C 

5  5.0 

7.C 

19.0 

8.0 

2.0 

7.0 

2.0 

2.0 

8.0 

2.? 

21.0 

6.0 

C  .  0 

*.r 

3  .  0 

2.0 

29.  0 

P.' 

27.0 

7.0 

0.0 

10. c 

3.C 

2.0 

50. r 

6.r 

17.0 

2.0 

0.    ) 

7.o 

1." 

2.n 

5  5.0 

5.r 

24.0. 

8.0 

2.0 

10.0 

2.0 

5.0 

8  5.  o 

5.0 

19.0 

4.C 

0.0 

6.0 

2.C 

3.C 

39.0 

6.0 

29.0 

0.0 

3.0 

13.0 

3.C 

6.0 

86.0 

-■  •  '. 

26.0 

3."^ 

1.0 

11.* 

2.0 

^.o 

95.0 

8.0 

4  3.0 

2.0 

1.0 

14.0 

1.0 

l.C 

104.0 

FFSLLTS 


Bl 


ESTIMATE  CF  P 


8. 1^782 
14.23561 

6.89968 
o.36C23 
5.3691c 
2. ^796^ 
2.'.  74 PA 


-0 


0. 

0. 

■0. 


B2 

0  2  389 
17268 
C.462^ 
01056 
C69P2 
0J627 
01857 


ESTIMATE  CF  CCVARIANCE  "ATRIX  SIG^fi 


4. 8^1 
7.131 
C.Q6C 
■C.6f C 
1.667 
r„zcc 


7.131 
30.781 
1.478 
o.  139 
6.593 
0.216 


i..  96  0 
1.^78 
7.538 
r.15o 
-0.312 
0.297 


-0.60  0 
139 
15^ 

1.  129 

697 
197 


n 


1.667 
6.5°3 
-'^.312 
C.697 
2.4Q1 
0.c26 


-l.o^?   -5.292   -0.841 


1.0  36   -0.06  2 


499 
216 
297 
197 

526 
°.5  7C 
0.602 


-.  . 
0. 

0. 


-1.942 

-5.2Q2 
841 
r36 
06  2 
6'^  2 
82< 


-o. 
1. 

-0. 

r 
v  . 

2. 


T  TEST  OF  THE  HYPOTHESIS  THAT  A 
CrMDQNENT  CF  THE  B2  VECTOR  EQUALS 
ZE&C  AT  A  ,1C  SIGNIFICANCE  LEVFL 


TFST  STATISTIC 

-1.C085C 

2.°054A 

-1.57087 
C. 02769 
^.12975 

-C. 77488 
l.C 3131 


RESULT 

ACCEPT 
REJECT 
ACCENT 
ACCEPT 
REJECT 
ACCEPT 
ACCEPT 


62 


DATA 


Y2 


ENGINE    MJMRER:     51675C 


Y3 


Yl 

5.0 

6.r 
e.o 

3.L 

«.£ 

5.C 

1C.  < 

8.0 


MJMBFP    TF    TATA     PCINTS    = 


Y4 


Y5 


Y6 


Y7 


XI 


1 1 .  0 

3.0 

0.0 

12. r 

2.r. 

2.0 

21.0 

3.0 

2.0 

0.0 

12.  r 

1.0 

2.0 

21.0 

13.0 

6.0 

1.0 

13.0 

3.0 

3.0 

27.0 

23.0 

7.0 

1.0 

12.0 

1.0 

5.0 

61.0 

13.0 

2.0 

o.O 

4.C 

r.c 

1.0 

86. ^ 

16.0 

6.C 

1.3 

10.0 

1.0 

4.0 

13.0 

15.0 

2.0 

1.0 

10.0 

c.c 

0.0 

lf.r 

34.  0 

8.0 

2.3 

12.0 

l.C 

3.r 

5^.0 

17.0 

6.0 

CO 

12.0 

2.G 

3.0 

22.0 

RESULTS 


«1 


ESTIMATE    CF    P 


5.^2833 

12.c*llc 

4.3579* 

r#cc2c^ 

13.C 3C2A 

1.74137 
2.^2961 


P2 

0,00514 
0.13476 
0.0^868 
0.0  32-"  7 
-0.063  3  5 
-0.01460 
0.00354 


FSTIMATE  CF  CCVARIANCE 


5.536 
15.*«3 
4.301 
0.897 
4.70C 
C.730 
1.764 


15. 443 

64.65^ 

16.711 

4.937 

11.656 

0.2^3 

6.840 


4.30  1 
16.711 
6.516 
1.272 
4.02  3 
1.188 
3.358 


CE     MATRIX     SIG*"\ 

0.897 

A. 700 

0.738 

1.764 

4.937 

11.656 

0.208 

6.84^ 

1.272 

4.023 

1.188 

3.358 

0.5  68 

0.716 

-0.025 

%376 

0.716 

5.559 

1.384 

2.324 

0.025 

1.384 

C.Q23 

0.736 

0.376 

2.324 

0.736 

2.594 

T  TEST  CF  THE  FYPCTHFSIS  THAT  A 
COMPONENT  CF  THE  82  VECTOR  ECUALS 
7FRC  AT  A  .10  SIGNIFICANCE  LEVEL 


TEST  STATISTIC 

f  .  1567C 
1.2022  5 
r. 24397 
f. 19728 
-1.92719 
-l.C Q022 
C.  15775 


RESULT 

ACCEPT 
ACCEPT 
ACCEPT 
ACCEPT 
REJECT 
ACCEPT 
ACCEPT 
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DATA: 


Yl 


Y2 


ENGINE    NtJMPER:     516575 


Y3 


Y* 


Y5 


Y6 


NUMPFP  OF  fATA  PCINTS  = 


Y7 


XI 


7., 

36.  C 

n    ■> 
<-..-> 

8.0 

l.: 

6.C 

22.'' 

5.' 

37.0 

0.0 

0.) 

12.'^ 

2.C 

2.0 

6*.p 

6." 

36.  ) 

3.C 

0.0 

11." 

2.C 

2.0 

5  8.r. 

17.  ) 

29.0 

4.0 

0.  ) 

9.0 

2.^ 

1.0 

28.0 

16.' 

3  5.0 

8.0 

1.0 

8.0 

4.0 

36. : 

1 6  .  ' 

5  J  .  0 

6.  3 

2.0 

17.  r 

4.0 

4.C 

112.^ 

10.* 

3  9.0 

fl.tl 

1.  ) 

10.  0 

4.0 

4.0 

85.0 

I  3  •  0 

34.0 

2.0 

0.  ) 

9.p 

5.C 

3.0 

63.^ 

C       I 

-  • 

25.0 

7.C 

1.0 

5.r 

o.c 

o.c 

27.0 

RFSULTS: 


HI 


ESTIMATE  CF  P. 


8.78.76C 

24.76978 

3.14465 

-r:. 16926 

5. .1866 

1. 15(  op 

2.18*2? 


82 

0.03202 

^.19532 
0.01951 
0.01313 
0.09021 

0.^3127 
?.C1276 


ESTIMATE  CF  CCVARIANCE  MATRIX  SIGVA 


27.619 

3.22") 

6.475 

1.022 

1.557 

5.329 

1.510 

3.220 

16.5  14 

-3.422 

0.356 

3.89e 

0.785 

6.478 

6.475 

-3.422 

11.253 

1.717 

-2.551 

1.427 

-1.0  85 

1.022 

0.356 

1.717 

%423 

-0.093 

0.016 

0.04? 

1.557 

3.898 

-2.551 

-•''.09  3 

2.936 

-0.3C1 

C.22  7 

5.828 

0.785 

1.427 

0.016 

-0.80  1 

2,630 

1.123 

1.510 

S.478 

-1.085 

n.r47 

0.227 

1.123 

3.671 

T  TEST  CF  THE  HYPOTHESIS  THAT  A 
CC^PONPNT  CF  THE  82  VECTOR  EQUALS 
ZFPC  AT  A  .10  SIGNIFICANCE  LEVEL 


TEST  STATISTIC 

C. 51881 

4.1  :93^> 
r. 49716 
1.72416 
'.49939 

1.64798 
0.55922 


RESULT 

ACCCPT 
REJECT 
ACCFPT 
ACCEPT 
REJECT 
ACCEPT 
ACCEPT 
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FATA 


Yl 


Y2 


ENGINF  NUMBER:  516534 


Y3 


Y4 


Y5 


Y6 


MJMBFP    OF    FATA    FCINTS    = 


Y7 


XI 


11. c 

27.0 

5.0 

0.0 

8.0 

1.0 

a.': 

17.^ 

14.U 

32.  J 

7.0 

1.0 

14.0 

1.0 

3.0 

19." 

15.P 

4  3.0 

6.0 

o.-) 

17.  0 

2.0 

4.r 

8C.^ 

15.0 

42.0 

7.0 

2.0 

15.0 

3.C 

i.C 

84.  n 

14.  n 

28.0 

7.0 

3.0 

8.0 

2.." 

3.0 

24.0 

18.< 

2  7.0 

1.0 

3.0 

7.0 

2.0 

2.0 

24.0 

14.0 

34.0 

10.0 

1.0 

9.0 

2.0 

6.0 

C.r 

17.0 

3  9.0 

1.0 

n.o 

15.0 

2.0 

5.0 

5^.0 

1*.  1 

38.0 

8.0 

0.  3 

10.0 

2.0 

3.0 

32. ^ 

FFSULTS 


Rl 


ESTIMATE    CF     E 


13.79272 
27.26^35 
6.4149S 
1.36?»5* 
7.34557 
1.36ri5 
4. 134P7 


B2 

0.02293 

0.18827 
-0.0167? 

-0.006  6  2 
0 

0 

-0.01812 


.10755 

.013  87 


ESTIMATE  CF  CCVARIANCc  MATKIX  SIGM4 


4* 

-0. 

-?. 

1. 

-r. 

c. 

-r. 


1C5 

638 
pep 

182 
425 

365 
298 


-0.638 
14.60  3 
6.490 
-3.101 
5.C73 
^.176 
3.199 


-3.898 
6.490 
If .260 
-C.495 
1.293 
0.  317 
0.430 


1. 
-3. 
-0. 


182 

101 
495 
802 
718 
0.383 
-1.313 


1. 

i. * 


-0.42  5 

5.07  8 

1.293 
-1.718 

6.C61 
-0.688 

1.617 


0.385 
^.176 
0.317 
0.383 
•0.668 
0.2*2 
•0,142 


-C.29  8 
3.199 
C.43C 

-1.31? 
1.617 

-0.142 
2.312 


T  TEST  OF  THF  HYPOTHESIS  THAT  A 
COMPONENT  OF  THF  82  VECTOR  ECUALS 
7ERC  AT  A  ,1C  SIGNIFICANCE  LEVEL 


TFST  STATISTIC 

C. 89167 
3.88142 
-f . 41124 
-C. 38869 
3.^4178 
2.22170 

-C. 93863 


RESULT 

ACCEPT 
OBJECT 
ACCEPT 
ACCEPT 
REJECT 
REJECT 
ACCEPT 
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DATA 


Yl 


Y2 


ENGINE  NUMBER:  509081 


Y3 


Y4 


Y5 


Y6 


NUVREP    OF    DATA    FCINTS    =    in 


Y7 


XI 


4. 

17.0 

6.0 

1.1 

5,0 

1.0 

1.0 

57.n 

l°.f 

2  9.0 

1.0 

0.0 

l.C 

3.C 

1.0 

91.0 

2.  ? 

19.0 

0.0 

1.0 

4.° 

0.0 

0.0 

63.0 

7.*'' 

32.0 

5.0 

0.0 

9.0 

2.0 

0.0 

18c.o 

e.{ 

3  9.  0 

Q.r 

0.  1 

9.0 

3.0 

2.^ 

209.0 

r     f 

-  .  - 

21.0 

0.0 

0.0 

5.."' 

0.0 

1.0 

237.0 

7.0 

27.0 

4.0 

CO 

5.^ 

3.0 

1.0 

148.0 

3.0 

26.0 

6.0 

0.0 

7.r 

l.C 

l.C 

110.0 

6.<. 

21.0 

8.0 

0.0 

6.0 

1.0 

0.0 

118.0 

9.0 

24.0 

2.0 

0.^> 

6." 

2.^ 

2.0 

106.^ 

*E SILTS: 


PI 


ESTIMATE    CF    R 


7.^20^2 
17.745<  * 
3.24261 
f.775Ff 
3.1  9P2T 
1.23716 
C.R56C C 


B2 

-0.00015 
0.:5879 
0.0C65S 

-C.00^37 
0.01973 
0.00275 
0.00261 


ESTIMATE  CF  CCVARIANCE  MATRIX  SIGVA 


25.50C 

15.287 

-3.121 

-1.1^03 

-5.738 

A. 377 

1.127 

15.287 

35.331 

7.989 

-0.818 

3.663 

5.959 

1.431 

-3.121 

7.c89 

11.689 

-0.158 

5.135 

1.226 

0.193 

-1.003 

-0.818 

-0.158 

0.  122 

0.055 

-0.226 

-0.r>5  3 

-5.738 

3.663 

5.135 

r.C5  5 

4.661 

0,00  2 

0.0*1 

^.377 

5.959 

1.226 

-'>.226 

0.00  2 

1.519 

■°.42n 

1.127 

1.431 

0.193 

-0.^5  3 

o.ou 

0.4  20 

0.584 

T  TEST  OF  THE  HYPOTHESIS  THAT  A 
COMPONENT  CF  THE  P2  VECTCR  EQUALS 
ZFRC  &T  A  .10  SIGNIFICANCE  LEVEL 


TEST  STATISTIC 

-r.C  ^5*6 
1.79515 
0.34506 

-2.27294 
1.65826 
0.^)507 
0.61910 


RCSULT 

ACCEPT 
ACCEPT 
ACCEPT 
REJECT 
ACCEPT 
ACCEPT 
ACCEPT 
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DATA: 


Yl 


Y2 


9.0 

21.0 

3  .  • 

11.0 

5.1 

2  2.0 

5.< 

19.0 

6.^ 

17.0 

9.0 

18.0 

7.0 

15.0 

6.( 

16.0 

6.0 

13.0 

ENGINF  NUMRER:  516592 


Y3 


O  .  v. 

5.C. 
^.^ 

5.0 

•  V 

/-.  n 

V  .   v.' 

6.0 
2.0 
4.C 


NUMBER  OF  r^TA  PCINTS  = 


Y4 


Y5 


Y6 


0.0 

12.0 

2.0 

1.0 

6.0 

2.Z 

0.0 

10.0 

A.Q 

1.0 

9." 

3.0 

1.0 

10.0 

2.0 

0.0 

11.0 

3.0 

l.J 

8.0 

2.0 

1.0 

10.  f 

2.C 

2.0 

4.0 

1.0 

Y7 


XI 


4.C 

89.0 

2.C 

30.0 

3.0 

59.0 

2.0 

83.^ 

4.C 

9  5.0 

5.C 

63.0 

4.0 

37.0 

0.0 

3  3.  r 

0.0 

34.0 

BFSU.TS 


01 


ESTIMATE  OF  B 


'-.6  66  21 
3. Of  392 
5.39C8C 
1.43G99 
5.2227f 
1.73505 
0.43833 


B2 

0.02678 
0.07297 
-0.0  3c2  3 
-0.01124 
0.06309 
0.01030 
0  .  C  3  8  3  5 


ESTIMATE  CF  CCVARIANCE  MATRIX  SIGMA 


3.66  8 

1.792 

-2.649 

-0.^18 

1.87C 

-0,3^8 

1.30  2 

1.7  =  2 

8.204 

-0.99  5 

-0.525 

1.222 

^.944 

-0.828 

-2. 649 

-0.995 

4.367 

0.548 

-2.358 

r.264 

-0.362 

-r.4ie 

-0.525 

^.5*8 

0.410 

-1.055 

-0.387 

-C.619 

1.67C 

1.222 

-2.35R 

-1,055 

4.196 

0.832 

1.227 

-G.3C8 

0.944 

0.264 

-0.3*7 

0.832 

0.775 

0.40c 

1.3C2 

-0.828 

-0.362 

-0.619 

1.227 

0.40  9 

2.579 

T  TEST  CF  THE  HYPOTHESIS  THAT  A 
COMPONENT  CF  THE  P2  VECTCR  EQUALS 
ZERO  AT  A  .10  SIGNIFICANCE  LEVEL 


TEST  STATISTIC 

l.( 2799 
1.P7332 
-1.3  8041 
-1.29033 
2.26479 
r. 65930 
!.7559c 


RESULT 

ACCEPT 
ACCEPT 
ACCEPT 
ACCEPT 
REJECT 
ACCEPT 
ACCEPT 
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DATA 


Yl 


Y2 


ENGINE  NUMBER:  515931 


Y3 


Y^ 


Y5 


Y6 


NUVBFP  OF  CATA  POINTS  =  10 


Y7 


XI 


18. o 

35.0 

0.0 

1.0 

11.0 

2.0 

3.0 

60.i. 

K,. 

27. r 

O.C 

0.0 

10.0 

3.0 

2.0 

36.0 

6.0 

34.  3 

3.0 

0.0 

11.0 

2.0 

3.0 

31.0 

ICC 

37.0 

7.0 

1.? 

14.0 

■     2.^ 

3.G 

120.0 

e.v 

^  ^.n 

5.0 

0.0 

10.0 

2.0 

2.0 

89.0 

12. c 

34.0 

6.G 

1.3 

ll.o 

0.0 

5.0 

89.^ 

5.C 

2  3.0 

5.0 

0.0 

8.0 

CO 

1.0 

1«.0 

10.  C 

33.0 

2.C 

l.n 

1 1 .  G 

CO 

1.0 

79.0 

6.C 

27.0 

c ) 

0.  1 

8."1 

0.0 

2.0 

51.0 

5.C 

21.0 

4.0 

o.o 

7.-> 

CO 

l.r 

9.0 

FFSUTS: 


FSTIMATE    CF    P 
Rl  B2 


6.-3751 
23.34^7 

i.cG7e: 

-C.18A?C 
7.48260 
C756-;e 
1.2°F6F 


C 04747 

j,  11599 
C. 02908 
0. 01004 

0.n4^96 
0  .  0  0  5  c  I 
O.C  172^ 


13.719 
8  .  704 

C966 
2.165 
1.698 
1.8  79 


ESTIMATE    CF    CCVAPIANCF    MATRIX    SIGM 


8.704 
13.569 
-2.737 
0.65  5 
3.884 
1.755 
2.225 


904 

.737 


-4 

_  n 

6.  ^fl« 

0.144 

•'..15- 

.  «96 

.2J9 


1.Q.5  6 
0.655 

-C  \<+H 

155 
.178 
-0.135 

10  2 


'  • 


2.165 
3.884 
-C15r 
0.178 
1.714 
0.982 
C47  3 


1.696 
1.755 
-0.8C6 
-0.135 
CQR2 
1.5£2 
1.192 


1.979 
2.225 
0.2GS 
C  1 0  2 

C  <>  7  F 
Cl°2 
1.339 


T  T£ST  CF  THE  HYPOTHESIS  THAT  A 
COMPONENT  CF  THE  P2  VGCTCP  ECUALS 
Z^RC  AT  A  ,1~  SIGNIFICANCE  IF.VFL 


TEST  STATISTIC 

1.37266 
3. 372^3 
1.22274 
2.72  90  5 
3.67812 
f. 50639 
1.59321 


RESULT 

ACCEPT 
PEJECT 

arcFPT 

REJECT 

REJECT 
ACCEPT 
ACCEPT 
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ENGINE  NUMBER:  5156r)P 


TATA: 


Yl 

£.( 
P.' 
a.  1 
7.'' 
LC.,0 
2.' 
*•< 

6. 


Y2 

39.0 
4  )./) 
31.0 
31.0 
37.0 
30.0 
21.0 
2.7  .,0 


Y3 

4.C 
1.0 

1  c .  0 

7.r 
l.«. 
3.C 

3.0 

4.0 

9.0 


Y^ 


Y5 


Y6 


Y7 


XI 


1.0 

16.0 

l.C 

3.-C 

2  2." 

0.0 

22. A 

2." 

4.0 

66.0 

0.0 

18.0 

l.C 

4.0 

58. r 

0.  0 

14.0 

1.0 

2.~ 

2P.0 

1.0 

21.0 

2.: 

6.C 

112.  C 

1.0 

15.0 

2.'^ 

3.0 

8  5.  0 

2.0 

12. C 

l.C 

0.0 

63.- 

0.  1 

12.  C 

3.0 

5.0 

36.^ 

1.  ) 

12. C 

0 . 0 

l.C 

27.0 

NU^PfcR  OF  DATA  PCIMS  = 


RESULTS: 


01 


FSTIMflTF    CP    P 


6.2C562 
3f  .  1  3 ^  7 C 

8.<  ^47 

T. 35927 
11.71p27 

0.785C? 

l.r5565 


R2 

.000  61 
.  •'•27  74 
.C611G 
0  .,00557 
0.07357 
0.01192 
".02817 


0 , 


.r* 


ESTIMATE  CF  CCVARIANCE  "ATRIX  SIGVS 


7.079 

1C.69* 

-3.229 

-0.473 

^.967 

C.4r2 

2.70C 


1 0  . 
39. 
-4. 
-2. 
18. 

0. 

6. 


684 
1Q7 

37h 

C18 
775 
274 
375 


_3 
-4 


—  r 

-i 


229 
3  74 
390 
359 
LI  9 


-1.335 

-r.585 


-.%^7  3 
-2.C 18 

-  ),35<3 

0.539 

-1.380 

-1.307 

-0.O7  3 


4.967 
13.775 
-1.119 
-1.330 
11.145 
0.^6° 
3.155 


0.452 

0.274 

-1.335 

-0.3F7 


^.741 
l.rt15 


2.7 
6.  37^ 

-c.ses 

-r.973 
3.  155 
1.015 
3.299 


T  TEST  OF  THF  HYPOTHESIS  THAT  A 
COMPONENT  Cc  THE  P2  VECTOR  ECUALS 
ZERC  AT  A  .10  SIGNIFICANCE  LEVEL 


TEST  ST4TISTIC 

-C.C1954 
^.37873 

-1..8C28C 
r. 64  801 
1.P3347 
1.  18436 
1.32551 


RESULT 

ACCEPT 
ACCEPT 
ACCEPT 
ACCEPT 
ACCEPT 
ACCEPT 
ACCEPT 
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ENGINE  NUMBER:  516380- 


CATA 


Yl 


Y2 


Y3 


Y4 


Y5 


Y6 


NUPBFP  OF  TATA  PCINTS  =  11 


FFSULTS: 


FSTIMATE  CF  P 
Bl  P2 


Y7 


T  TEST  CF  THE  HYPOTHESIS  THAT  A 
COMPONENT  CF  THE  B2  VECTOR  ECUALS 
ZERC  AT  $    .10  SIGNIFICANCE  LEVEL 


XI 


5.0 

19.0 

6.C 

0.0 

9.0 

1.0 

4.0 

6^.° 

11." 

3  5.0 

3.n 

1.0 

15.0 

l.C 

1.0 

5  8.0 

18.  0 

43.  Q 

5.0 

2.0 

13.  n 

2,0 

5.0 

5  2.0 

A.r 

?4.0 

7.0 

2.0 

10.0 

l.C 

l.C 

12.0 

6.0 

27.0 

9.C 

l.o 

10.0 

1.0 

l.C 

2<;.C 

6.( 

25.0 

1.0 

1.  3 

12.  0 

1.0 

l.C 

41.0 

5.< 

26.0 

6.G 

0.  ) 

12.0 

1.0 

2.0 

66.° 

c     r. 

-  . 

19.0 

7.0 

3.0 

8.0 

1.0 

5.0 

6.0 

*••   ) 

16. G 

5  .  0 

n.  3 

4.^ 

1.0 

0.0 

39.0 

K.^ 

1  5  .  0 

5.0 

1.0 

5.C 

1.0 

o.c 

31.0 

3. 

19.0 

2.A 

0.0 

6.0 

1.0 

2.C 

18.0 

3.26T27 

0. 

08528 

17.74008 

0. 

17684 

5.94440 

_  o 

02279 

1.9773'' 

• 

C.26^9 

6.2420C 

0. 

08577 

C.958C6 

0. 

00  3^5 

1.P7213 

0. 

G  ^341 

FSTIMATE 

CF    COVARIANCE    MATRIX    SIGMA 

16.545 

29.482 

CO  57 

2.681 

6.526 

1.145 

3.756 

28.482 

65.145 

0.463 

4.880 

2C565 

i.7e5 

5.392 

O.f  57 

0.463 

5.864 

o.5'7 

-0.049 

0.027 

0.924 

2.681 

4,88  0 

0.5T7 

0.801 

1.909 

0.153 

0.92  9 

6.526 

20.665 

-C.C4  9 

1.909 

1 0 .062 

1.255 

1.644 

1.  1^5 

1.785 

0.02  7 

0.153 

C  2  5  5 

0.095 

C.  328 

3.7f6 

5.392 

C.924 

A.C29 

1.644 

C.328 

3.772 

TEST  STATISTIC 

1.34267 
1.40306 
-f. 60264 
-1.8o709 
1.73153 
C. 73587 
r. 11256. 


RESULT 

ACCEPT 
ACCEPT 
ACCEPT 
REJECT 
ACCEPT 
ACCEPT 
ACCEPT 
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TATA 


VI 


Y2 


IF... 

41.0 

IP.' 

4  3.') 

17-.' 

*■  5. 'J 

12,0 

5  6,  0 

12.' 

37.  ) 

U  .  • 

3?.'1 

C." 

34,0 

1  Pa- 

*-■ *> .  ') 

ll. 

3  6.': 

ENGINE  NUMPER:  515515 


Y3 


NUN'RcP  OF  TATA  POINTS  = 


YA 


Y5 


7.0 

C.  1 

13.; 

1.0 

1.0 

11.0 

3.0 

C  .  0 

11.0 

9.~ 

1,0 

12.  r 

7.0 

O.'l 

10.C 

3.'" 

1.0 

8  .  0 

1.0 

1.  J 

13.^ 

2.1 

0.0 

15. C 

6.o 

0.0 

a .  o 

Y6 

2.0 

3.0 
2.0 
2.C 
2.C 

2.0 
l.C 

2.0 

2.0 


Y7 


XI 


6.C 

42.: 

5.0 

55.0 

'-.C 

5  8.0 

4.C 

37.0 

l.C 

24.0 

0.0 

24.0 

3.0 

62.- 

?..r 

76.^ 

2.^ 

3.0 

RESULTS: 


Qi 


ESTIMATE  CF  P 


If.  12277 
32. 5r  7£* 

7./ -C  92f 

r . 175CC 

7.577F1 
2. .70oq 
0.«=7r21 


P2 


n 

/  a 


0  8896 
.16124 


0.C7739 
0,0 009 6 
C.C86^S 
3*0016  8 
".055R? 


FSTIVATR  CF  CCVAPMNC-  MATRIX  SIG^A 


11.693 
14,C£6 
2.  153 
-l.*73 
C.46c 
1.375 
3.9^4 


14.r66 

20.499 

1.440 

-2.194 

-.057 
1.733 
4.^62 


2.153 

-1.473 

*.465 

1.375 

3.944 

1.440 

-2.194 

:.r57 

1.733 

4.^6  2 

6.^29 

-n.044 

2.5Q1 

-0.^77 

2.147 

0.044 

<^.2P5 

-~'.^46 

-0.144 

-C.397 

2.591 

-0.rA6 

1.803 

-0.2 CO 

1.182 

0.077 

-'-.144 

-0.200 

0.2  84 

^.342 

2.1^7 

-0. 397 

1.182 

n.342 

2.713 

T  TEST  CF  THE  HYPOTHESIS  THAT  A 
COMPONENT  CF  THE  62  VECTCR  ECUALS 
ZERC  AT  A  .10  SIGNIFICANCE  LEVEL 


TEST  STATISTIC 

1.6  3  ISC 
2.30081 

-l.c717? 

-0. 11594 

^.14212 
-f. 23330 

2.13941 


RESULT 

ACCEPT 
RE  IECT 
REJECT 
ACCEPT 
REJECT 
ACCEPT 
PEJECT 


71 


CATA 


Yl 


Y2 


e.r 

16.  a 

11.  t 

29.') 

13.  C 

21.0 

bmV 

20.0 

20.0 

10.' 

22.0 

13. C 

24.0 

5.0 

23.0 

5  .  0 

15.0 

ENGINE    NIJMRER:     516731 


Y3 

6.o 

6.0 
fl.C 

6.0 
8.0 
5.0 
5.0 

7.0 


NUMBFP  HF  DATA  PCINTS  = 


Y4 


Y5 


Y6 


1.^ 

14.0 

2.0 

2.0 

16.0 

3.0 

2.0 

9.0 

1.0 

3.? 

9.r> 

0.0 

o.o 

8.0 

1.0 

1.0 

13.0 

2.0 

1.0 

7.0 

1.0 

2.0 

9.0 

2.e 

1.0 

8.0 

1.0 

Y7 


XI 


5.0 

5^.1. 

6.C 

8  3.  r 

3.0 

57.0 

5.C 

33.0 

1.0 

87.0 

3.^ 

27.0 

0.0 

2  2.  C 

4.0 

99.  0 

0.0 

31. r 

RESULTS: 


SI 


FSTIMATE  CF  P 


10.3965' 
18.453"? 
4.tC35r 
1.45355 
c.  1PP3- 
0. 67471 
1.45051 


P2 

C3564 

C2G55 
C0017 

C209c 

"l  405 
0.02829 


■*>  . 

-0. 

0. 


ESTIMATE  CF  CCVAPMNCE  MATRIX  SIGMA 


11.962 
9.  673 
G.45C 

C.312 
3.521 
1.224 

1.248 


9.B73 
17.879 
-1.115 
1.231 
3.984 
1.429 


0.450 

-1.115 

0.817 

0.021 

-1.399 

-0.523 


2.261       -L.369 


0.312 

1.230 

0.021 

0.889 

0.  384 

-0. 109 

1.433 


3.521 
3.984 

-1.3°9 
<">.384 

11.011 
2.243 
5.^36 


1.224 
1.429 
-'■,.523 
-0.109 
2.243 
0.700 
0.621 


1.24P 
2.261 
-0.369 
1.433 
5.436 
0.621 
4.951 


T  TEST  CF  THF  hYPCTHESIS  THAT  A 
COMPONENT  CF  THE  82  VECTCR  ECUALS 
Z^RO  AT  A  .10  SIGNIFICANCE  LEVEL 


TEST  STATISTIC 

-0.84223 
C.c3805 
2.67191 

-C. 01442 
r. 51491 
1.37245 
1.(392^ 


RESULT 

&CCEPT 
ACCEPT 
REJECT 
ACCEPT 
ACCEPT 
ACCEPT 
ACCEPT 
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ENGINE    NIJMRER:     515447 


TATA 


vi 

1 G  .t 

1?.. 
P. 
7. 
P.- 

r     r, 
.  # 

A. 

6.  : 
12. ° 

12. i 


Y2 


25. 
2". 
26. 
27.0 

21.0 

27. 


'j 


39.0 


Y3 

9.  • 
2." 
4.  ■ . 
4.C 
1.0 
3.0 
4.0 
7,~ 
3.r 
2.C 


Y^ 


Y5 


• 

11." 

3.° 

15.' 

o.  "> 

9.'" 

1.0 

11." 

?.  ■) 

10..': 

l.  i 

<?.c 

1.0 

•3." 

1.0 

c>.<" 

12.  ? 

9  .  "* 

Y6 

l.( 
2.' 

i.C 

l.r 
l.r 

1 . 0 
l.r 


4 .  0 


Y7 


XI 


2.r' 

6  5.  n 

2.^ 

82.^ 

0.0 

46." 

?.r 

54,  " 

3.^ 

6  v .  c 

o.c 

71.0 

CO 

2  5. '-1 

3.0 

26,0 

2.C 

82.0 

".r 

7?"     •" 

NU>'QFP     OF    TATA     PCI  MS     =     lf 


PE  SULTS 


PI 


t ST! MATT    CF    P 


1.992C1 

16.25591 

6.f^£2fc 

-r.t  C2C5 
t.',2cS2 

c.  P4f  94 

1.Z562P 


P2 

110 10 
222^1 

C4710 

;13P0 

'^•64  8? 
J25."7 

C02^7 


estimate;  cf  ccvariance  matrix  SIGMA 


4,4(4 

3.621 
C.  P'  5 

-0     3  c  7 

dp2  5 
2.0  65 

C673 


3.621 
2  1 .  9 1  4 
4,162 
1.496 
7.355 
0.  2C3 
3.519 


0.845 
'•.162 


-r 


•  57o 

.34  7 

J.P2^ 

--C7P7 

C729 


-0.357 

1.4C6 

-•  =  .847 
1,111 
1.133 

-1.^86 
a.  7o<3 


C  .  8  2  5 
7.358 
0,825 
1.133 
3,242 
-',-.83  2 
1.4  5C 


2.065 
C2C3 
-C787 
-1.086 
-C832 
3,97G 
-1.054 


C  6  7  3 
3.519 
0.729 
C7^9 
1.45C 
-1.P54 
1.797 


T  TEST  CF  THE  HYPOTHESIS  THAT  A 
CCVPGNENT  CF  THE  82  VECTCR  EQUALS 
ZERO  AT  A  .10  SIGNIFICANCE  LEVEL 


TEST  STATISTIC 

3.20222 
2.89091 
-1.21977 
C .7^880 
2.19847 
C. 76801 
r.11243 


RESULT 

REJFCT 

REJECT 
ACCFPT 
ACCEPT 
REJECT 
ACCEPT 
ACCEPT 
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ENGINE 


NUMBER:  516048 


TATA 


Yl 


Y2 


Y3 


Y4 


Y5 


Y6 


Y7 


MJ^PEP  PF  CATA  PCINTS  = 


XI 


6 .  '■,- 

24.0 

1.0 

0.3 

8.0 

3.0 

2.C 

30.^ 

G.c 

*6.0 

4.^ 

2.0 

12." 

2.0 

3.0 

120.0 

4.r 

24.0 

6.C 

1.0 

6.0 

0..G 

0.0 

90.  * 

c.  0 

46.0 

5.0 

1..) 

13.0 

2.0 

4.0 

96.0 

l.i 

19.0 

3.0 

2.  3 

5.0 

1.0 

1.0 

59.  r 

6.  " 

3  9.0 

6.C 

1.0 

ll.Q 

l.C 

1.0 

5  5.0 

4.o 

29.0 

3.0 

2.0 

9.: 

1.0 

3.0 

33.^ 

C.< 

•52  .0 

8.>3 

1.0 

16.  C 

1.0 

2.0 

96.  C 

2. 

29.0 

c.c 

0.? 

9  .  0 

1.0 

0.0 

2°.r' 

FESULTS: 

ES 

TIMATE    CF    B 

Bl 

P2 

1.94411 

u. 06004 

19.C5C2* 

0.22459 

r.!:3192 

C. 05134 

C.4833C 

0  .  c  c  9  2  9 

6.19133 

0.C53C 9 

f .2946? 

0.CK44 

C. 8364 7 

0. CI  393 

FS 

TIMATE    CF    CCVARIANCE     MATRIX    SIGMA 

4.123 

14.290          0.952 

-0.304 

5.385 

C.319 

1.9^1 

14.29C 

92.544          5.^2P 

-1.773 

30.5*6 

3.20  6 

6.096 

r.9c2 

5.428         3.967 

0  .  r-  8  8 

2.277       - 

-C.418 

-^.368 

-0.3C* 

-1.77  3         0.088 

0.585 

-C.616 

*.158 

ft.  290 

5.385 

30.5*6         2.277 

-0.616 

10.5  2C 

0,843 

2.251 

C.319 

3.206      -r.418 

0.158 

0.843 

0.42  8 

C.523 

1.9C1 

6.C<5  6      -0.368 

0.290 

2.251 

0.523 

1.967 

T  TEST  CF  THE  HYPOTHESIS  THAT  A 
COMPONENT  CF  THE  82  VECTCP  ECUALS 
ZERO  AT  A  ,1ft  SIGNIFICANCE  LEVFL 


TEST  STATISTIC 

2.E3524 
2.23854 
2.^7146 
1.16508 

1.56945 
l.c2994 
C.<=5259 


RESULT 

REJECT 
REJECT 
REJECT 
ACCE°T 
ACCEPT 
ACCEPT 
ACCEPT 
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ENGINF  NUMBER 


16*63 


CATA: 


Yl 


Y2 


Y3 


MJMRER  CF     rATA  PCINTS  = 


Y4 


Y5 


2. 

11.0 

3.C 

0.0 

6.0 

7.r 

20.0 

0,0 

0.0 

12.0 

6.r 

18.0 

2.0 

0.0 

10.0 

8*0 

19.0 

CO 

1.0 

CO 

2.0 

13.0 

8.0 

1.  ■) 

8.0 

S<.0 

23.0 

l.o 

0.0 

14.  n 

4.G 

17.0 

9.0 

u« 

9.0 

5*0 

2%0 

5.0 

1. ) 

9.0 

7.r 

32.0 

2.C 

0*0 

14.0 

Y6 


Y7 


XI 


6.0 

2.0 

CO 

1.0 

4.0 

71.0 

1.0 

2.*" 

34.0 

1.0 

4»n 

40.  n 

o,r 

2.0 

2  7.0 

1.0 

2.^ 

0**5 

2.0 

4.0 

2^.0 

7.0 

2.0 

17.0 

7.0 

4.0 

6  5.0 

FfSlLTS: 


1 


FSTIMATE  OF  P 


<,.  53252 
15.26°6* 
4.63345 
f . 5^999 
e. 66^85 
3.1106c 
l.=4C8(. 


R2 

0.03312 

0*12*72 
-0.0420  9 
-0.00309 

C. 04695 
-0.00718 

G.C3*4C 


ESTIMATE  CF  COVARIANCE  MATRIX  SIGMA 


6.387 

9.  302 

-5.951 

-0.3  87 

5. 089 

-1.321 

0.356 

9.30  2 

30.475 

-5.324 

-0.850 

12.190 

6.961 

C.450 

5.951 

-5.32^ 

11.299 

1.  145 

-3.342 

-0.027 

0.253 

C  .  3  8  7 

-0.P5O 

1.145 

0.311 

-0.673 

-0.238 

•C  1 3 1 

5.0  89 

12.19* 

-3.342 

-0.673 

6.829 

-0.170 

-O.OIO 

1.321 

6.961 

-0.027 

-0.238 

-0*170 

9,519 

-0.00  2 

C.  356 

0.4  50 

0,253 

0.131 

-o.oio 

-0.00  2 

n.*r>s 

T  TOST  CP  THE  HYPOTHESIS  THAT  A 
COMPONENT  CF  THE  82  VECTCR  ECUALS 
ZERC  AT  A  .10  SIGNIFICANCE  LEVFL 


TEST  STATISTIC 

C.c2937 
1.60218 

-0.83798 
-C. 39360 

1.27409 

-0.16505 

2.76992 


RESULT 

ACCENT 
ACCEPT 
ACCEPT 
ACCEPT 
ACCEPT 
ACCEPT 
REJECT 
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DATA: 


Yl 


Y2 


ENGINE  NUMBER:  539064  , 
Y3      Y4      Y5      Y6 


ESTIMATE  CF  P 


Bl 

7.72752 
25.42232 
6.27336 
f  .36*8' 
6. 1518? 
3.73C31 
1.86632 


B2 


0.02695 
o .  r  7  P  8  5 
0.0G297 
C .00576 
3. 06053 
■C. 00820 
O.C1630 


Y7 


XI 


It.? 

^5.3 

l.r 

1 . 0 

2  2.0 

2.C 

4.0 

240.0 

ll.r 

31.  1 

2.0 

1.  ) 

11.0 

3.0 

5.C 

10  3.0 

9..C 

3*3.0 

6.0 

0.0 

1 3 .  C 

4.0 

2.0 

90.0 

10.0 

3  5.0 

7.0 

1.0 

11.0 

4.0 

4.0 

67.  n 

1C.0 

24.3 

7.0 

l.n 

8.0 

5.C 

1.0 

30.0 

14.0 

^2.0 

P.r 

2.0 

19. r 

2.0 

6.0 

21G.0 

14. 

'3.0 

8.0 

3.0 

19.  C 

2.0 

6.0 

217.  0 

11." 

35.0 

9.0 

0.0 

13.  n 

2.0 

5.0 

150.0 

8.  3 

31.0 

5.C 

1.3 

10.  c 

0.0 

3.0 

60.0 

NUMBER  OF 

D£TA  PCINTS  = 

g 

RESULTS: 

FSTIMATt  OF  CCVARIANCE  MATRIX  SIGMA 


C.9CC 

-1.163 

C.290 

0.483 

-0.10  8 

0.8  86 

0.151 

1.163 

8.229 

0.152 

-0.098 

2.517 

-0.063 

-O.150 

C.29C 

0.152 

8.639 

0.561 

-1.0  64 

0.358 

1.043 

t.^63 

-0.09  8 

T.561 

^.75^ 

C  •  C  5  7 

-0.O59 

0.351 

c.ire 

2.517 

-1.064 

0.057 

1.626 

0.381 

-1.00  4 

r.686 

-O.06  3 

C.  35  8 

-0.059 

0  .381 

2.113 

-0.518 

C.151 

-".150 

1.043 

0.361 

-1.CP4 

-0.518 

1.619 

T  TEST  OF  THE  HYPOTHESIS  THAT  A 
COMPONENT  CF  THE  P2  VECTOR  ECUALS 
7FRT  AT  A  .10  SIGNIFICANCE  LEVEL 


TFST  STATISTIC 

5.91299 
6. (,7753 

-(.22:25 

1.44270 

10.36179 

-1.23196 

2.79686 


RESULT 

REJECT 
REJECT 
ACCEPT 
ACCEPT 
REJECT 
ACCEPT 
REJECT 
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DATA 


Yl 


1 


Y2 


12.  ? 

3  3.0 

9. 

*  vo 

ic; 

4  3*0 

9.' 

*?.  ) 

e.  • 

30.0 

7VC 

2  3.0 

ll.i 

42*0 

n:*o 

45.0 

'O. 


ENGINE  NUMBER:  ^24646 


Y3 


2. 

L> 

7. 

r 

1. 

r 

0. 

r> 

2. 

/ . 

s 

S 

w  • 

3. 

Q 

6. 

L. 

NUN'BER  np  DAT*  FCINTS  = 


Y4 


Y5 


r. ) 

1 1.0 

0-.G 

13.0 

2.1 

11.0 

0.0 

9.0 

0.0 

10.0 

1.0 

8  •  0 

0.0 

12.0 

0.0 

12.  r 

1P.O 

Y6 

2.0 
2.C 
2.0 
3.0 
l.C 
1.0 
2." 
2.o 
5  .  0 


Y7 


XI 


4.0 

3  2.0 

6.C 

32.  n 

5*0 

27.0 

2*0 

22.0 

1*0 

S..0 

1.0 

8.0 

5.^ 

63.^ 

3.C 

72. r 

7.0 

126.0 

F  0  S  L  L  T  S 


Bl 


EST  I'M  ATt  CF  P 


B2 


7.7630  4 
31*55843 
2.17532 
0.-55916 
6.36977 
£.,110359 
2.16722 


0  5  6  29 
19645 

CO  36  5 

0590 

562 

02581 

0,03947 


Q  . 

.-1 


V 


1.2*5 
1.345 
-0.71P 
C .  C  2  3 
0.377 
C.260 
C.6Q3 


ESTIMATE  CF  CCVARIANCE  MATRIX  S I G w  A 
5    1.343   -3*719    ^.C23    ".377 


1.34  j> 
13.20'J 
0.2^2 
0.146 
0.747 
1*91.0 
2.558 


-:.71Q 
0.202 
7.692 

-0,67  9 
0.546 

-1.107 
1.287 


,.  146 


-'^67  9 


.377 
.747 


%514 
0.0  20 

0*01.3 
0,3  3fc 


0.54  6 
0*020 
2.0L9 
U  .  3  1 2 

1.5  64 


0.  26C 
1.910 
-  1 .  1  C  7 
G.013 
°.31  2 
0.554 
0.244 


0.69  3 
2.558 
1.287 
0.336 
1.564 
0.244 
2.801 


T  TEST  CF  THE  HYPOTHESIS  THAT  A 
COMPONENT  CF  THE  62  VECTCR  ECUALS 
ZE*C  AT  a  .10  SIGNIFICANCE  LEVEL 


TEST    STATISTIC 

5.2^951 
5.P^30  0 
C  .  1 4 1 1 0 

-r. 8 8370 
5.71138 
3.72125 
2.53127 


RESULT 

REJECT 
REJECT 
ACCFPT 
ACCEPT 
REJECT 
REJECT 
REJECT 


77 


ENGINE  NUMBER:  516627 


DATA: 


Yl 


Y2 


Y3 


Y4 


Y5 


Y6 


MJM9EP  OF  TATA  FCINTS  =  11 


PEStLTS: 


PI 


ESTIMATE  CF  P 
62 


4.  95*-7<' 

0.0612*5 

1.82213 

0.18915 

?. ^9675 

?. 02352 

C.?r373 

0.00378 

6.62466 

0.090  17 

1.672C2 

0  •  C  0  9  9 1 

l.*173P 

0. 01765 

Y7 


FSTIMATE  CF  CCVARIANCE  MATRIX  SIC-MA 


T  TEST  OF  THE  HYPOTHESIS  THAT  A 
COMPONENT  CF  THE  62  VECTCR  ECUALS 
ZERC  AT  A  .1C  SIGNIFICANCE  LEVEL 


XI 


8.' 

32.0 

9.0 

0.0 

7.0 

3.0 

3.C 

24.  C 

*.  ) 

18.0 

1.0 

0.0 

7.C 

3.0 

3.0 

10.0 

11. C 

37.0 

9.0 

2.0 

14. C 

3.0 

4.0 

86.0 

9.C 

27.0 

2.0 

0.0 

9.0 

2.0 

4.0 

2*.0 

5.C 

2  3.0 

2.0 

1.0 

7.0 

1.0 

o.o 

3.0 

ICC 

32.0 

3.0 

0.0 

13.0 

2.r 

3.0 

63. C 

1 C  .  0 

39.0 

2.0 

0.0 

13.0 

2*0 

A.o 

68.0 

8.0 

3  5.0 

n.o 

0.0 

13.0 

2.0 

3.* 

68.  « 

6.i 

28.0 

2.° 

1.0 

11. C 

1.0 

l.r 

23.0 

f.r 

27.0 

2.0 

0.0 

11. r 

1.0 

l.o 

36.0 

6.'J 

30.0 

2.C 

o.o 

1  2  .  0 

3.0 

4.0 

60.  0 

2.  669 

3.279 

2.382 

".116 

-0.339 

-J. 021 

C.172 

3.  278 

11.578 

3  .  9 1  J 

0.029 

-0.452 

-0.571 

-2.038 

2.382 

3.910 

9.631 

^.995 

-1.871 

1.125 

1.27^ 

C  .  1 1 6 

0.02  9 

0.995 

1.493 

0.094 

-0.072 

-0.107 

0,339 

-0.452 

-1.871 

0.09* 

1.157 

-0.718 

-0.746 

0.021 

-0.571 

1.12  5 

-0.072 

-r;.718 

0.684 

0.922 

0.172 

-2.C88 

1.274 

-0. 107 

-0.746 

0.922 

2.463 

TEST  STATISTIC 

3.27932 
4. 85956 
C. 66246 
C. 47103 
7.32913 
1.04718 
C.9934r 


RESULT 

REJECT 
REJPCT 
ACCEPT 
ACCEPT 
REJECT 
ACCEPT 
ACCEPT 
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ENGINE    NIJMPFR:    516101 


TATA 


Yl 


Y2 


Y3 


Y4 


Y5 


16.  ." 

44.  3 

lfc.n 

40.  0 

p. 

17.0 

l  r  :■ .  r 

33.0 

1 1 . r 

36.0 

6,1 

?9.0 

If  .f 

37.0 

7.   > 

21.0 

6." 

24.0 

7.r 

22.0 

u.' 

?  5  .  0 

3.0 

0.0 

2  3  .  C 

0.0 

0.  ') 

23.0 

1.0 

3.0 

8.0 

9.C 

UO 

19.0 

3.0 

A       '> 

21.0 

9.0 

0.0 

17.' 

1.0 

C.J 

16.0 

5.0 

2.0 

in/i 

2.0 

0.0 

11.  r 

8.0 

2.0 

9.0 

7.0 

0.0 

1 5  .  0 

Y6 

3.° 
1.0 

n.c 

1.0 

1.0 
1 .  0 
2.0 
1.0 

i.o 

l.C 
1.0 


MjfBEP    nF    DATA    PCI  MS    =    11 


PESUl  TS: 


31 


ESTIMATE    TF    P. 

B2 


P.(  1412 

0.02141 

24.5572* 

•>.<:6576 

5.  14212 

-'.0097  3 

r. 71736 

0.00012 

11.69686 

0.04924 

1.276(3 

-0.00118 

1.73C59 

0.00223 

ESTIMATE    CF    CGVARIANCE    «AT3IX    SIG*>A 


Y7 


XI 


4.C 

8C:.C 

l.C 

12P.C 

4.0 

12  8.0 

l.C 

192.  r' 

2.0 

114.? 

'"• .  0 

56.0 

3  .  C 

90.  C 

l.C 

23.^ 

1.0 

23.0 

4.0 

11.0 

1.0 

35." 

13.487 

26.779 

-2.465 

-2.767 

16.654 

2.040 

0.132 

26. 77^ 

69.565 

-6.T69 

-7.645 

38.864 

5.651 

T.676 

-2.46c 

-6.069 

12.388 

n.348 

-0.676 

-0.344 

-1.771 

-2.767 

-7.645 

1.348 

1.35  3 

-4.921 

-0.494 

0.635 

16. 6  54 

33.E64 

-0.676 

-4.921 

25.366 

2.733 

-1.545 

2.L4C 

5.651 

-0.344 

-0.494 

2.733 

0.621 

0.36  3 

C.132 

0.676 

-1.771 

'^,6  3  5 

-1.545 

0.363 

2.750 

T  TEST  CF  THE  HYPOTHESIS  THAT  A 
CCVDONENT  CF  THE  B2  VECTCP  EQUALS 
7ERC  £T  A  .10  SIGNIFICANCE  LEVEL 


TEST    STATISTIC 

1.C4745 

1.41632 

-C.*9665 

( .( 1014 

1.7563* 

-r. 25839 

0. 24171 


RESULT 

ACCEPT 
ACCEPT 
ACCEPT 
ACCE°T 
ACCEPT 
ACCEPT 
ACCEPT 
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DATA 


Yl 


Y2 


7.  ? 

25.0 

16.' 

29.0 

c.  "* 

35.  o 

5.. 

31.  Q 

9. 

3  1 .  0 

ic. : 

2  9.0 

6.i. 

23.0 

7.. 

24.0 

7.0 

2  )•  J 

ENGINE  NUMBER:  516132 


Y3 

9.0 
3.0 
6.0 
7.C 

4.r 

7. : 

3.0 
2.0 

1 0  •  0 


I^UVRER  OF  CAT A  PC  I  NTS  = 


Y4 


Y5 


2.1 

B.r 

0.0 

9.0 

2.1 

11.  P 

1.0 

9.C 

0.0 

9.P 

C.o 

11.  r 

0.0 

7.C 

0.0 

7.P 

0.0 

5.r 

Y6 


Y7 


XI 


1.0 

2.r 

26.0 

2.C 

0..C 

100.1 

2  .  0 

4.r 

139.0 

Z.C 

2.r 

73.0 

2.0 

c.c 

107.0 

1.0 

o.c 

5P.n 

2.C 

2.0 

44.0 

1.0 

1.0 

58.^ 

1.0 

2.0 

35.  n 

RESULTS: 

ESTIMATE  CF  P 

PI  R2 


c.c^c^l 

3.04049 

9. 6876* 

0.10908 

7.67937 

-0.02830 

r.?53C= 

0.00425 

t.:243c 

r'.C34C3 

C. 87049 

0  .  0  C  c  6  3 

1.15398 

0.00408 

FS 

FIMATE    CF    CCVAPIANCE    MATRIX    SIGV 

A 

^.420 

-1.576       -1.412 

-1.020         j. 551 

-0.225 

-2.66C 

-1.576 

6.789          2.53° 

1.229          3.546 

5.008 

0.^36 

-1.412 

2.533          7.P66 

1.4  3P          C.86  8 

-C.327 

1.660 

-1.C2C 

1.229         1.430 

0.860         0.45  2 

-0.034 

0.041 

0.55] 

3.546         0.668 

''.452          2.472 

-C126 

-0.333 

_p 


-2.660 


22 


3.009 

0.036 


_  r 


327 
1.66  5 


-n 


034   -0.126    J.  170 
0.941   -0.333    Q.^48 


C04P 
2.^0  5 


T  TEST  CF  THE  HYPOTHESIS  THAT  A 
COMPONENT  CF  THE  e2  VECTCR  ECUALS 
ZFRC  AT  A  .10  SIGNIFICANCE  LEVEL 


TEST  STATISTIC 

1.39313 
4  .  l 2 1 2  7 
-1.06572 
C  8443 
2.28583 
2. A  70  5  3 
C. 30463 


RESULT 

ACCEPT 
REJECT 
ACCEPT 
ACCEPT 
REJECT 
REJECT 
ACCEPT 
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F.   CCMFARISCN  CF  ESTIMATED  PARAMETERS  AMONG  ENGINES 


ESTIMATEC 

PARAMETERS 

FCR  CHROMIUM 

SERIAL  # 

MEAN 

VARIANCE 

516726 

3.80C0 

7.7333 

515355 

4.2222 

10.6944 

516678 

5.0000 

7.2  500 

515C56 

4.5556 

7.2778 

5C9C90 

3.1818 

6.  1636 

524739 

2.5714 

13.3406 

5163CS 

3.2222 

12.4444 

516C94 

5.5556 

11.5278 

516S4S 

4.1000 

8.7667 

51675C 

4.6667 

5.75C0 

516575 

4.2222 

10.1944 

516534 

5.7778 

9.1944 

509C81 

4.1CCC 

10.5444 

516592 

3.1111 

4.8611 

515S31 

3.2000 

6.8444 

5156C8 

4.6667 

1C.75CC 

516380 

5.09C9 

5.49CS 

515515 

4.3333 

8.75C0 

516731 

6.2222 

1.4444 

515447 

3.9000 

5.8778 

516C46 

4.0000 

6.5C00 

516C63 

3.3333 

ll.OCCC 

509064 

5.8889 

7.6111 

524646 

2.33  3  3 

6.7500 

516627 

3.09C9 

9.09C9 

516101 

4.3636 

11.4545 

516132 

5.6667 

8.CCCG 

ESTIMATED 

PARAMETERS 

FCR  SILVER 

SERIAL  U 

MEAN 

VARIANCE 

516726 

1.60C0 

0.7111 

515355 

1.0000 

1.2500 

516676 

C.7778 

0.4444 

515C58 

0.5556 

0.2778 

5C9C9C 

0.2727 

0.2182 

524739 

1.5714 

2.4176 

516309 

0.4444 

0.5278 

516C94 

0.6667 

C.75C0 

516949 

l.OOCO 

1.1111 

51675C 

0.6667 

C.5CC0 

516575 

0.5556 

C.5278 

516534 

1.1111 

1.6111 

509081 

0.2000 

0.1778 

516592 

0.7778 

0.4444 

5  15931 

C.40CC 

0.2667 

515608 

0.6667 

0.5000 

51638C 

l.OOCO 

l.OCCO 

515515 

0.3333 

C.25C0 

516731 

1.4444 

0.7778 

515447 

C.8CC0 

1.0667 

516048 

1.1111 

0.6111 

516C63 

C.4444 

0.2778 

509C64 

1.1111 

C.86U 

524646 

0.3333 

C.5CC0 

516627 

C.3636 

0.4545 

516101 

C.7273 

1.2182 

516132 

0.5556 

0.7778 
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ESTIMATED  PARAMETERS  FOR  MAGNES. 


SERIAL  * 


MEAN 


VARIANCE 


516726 
515355 
516678 
515058 
509C9C 
524739 
516309 
516C94 
516949 
516750 
516575 
516534 
509081 
516502 
515931 
515608 
516380 
515515 
516731 
515447 
516048 
516063 
509C6  4 
524646 
516627 
516101 
516132 


1.600C 
1.3333 
2.2222 
0.7778 
1.1818 
2.5714 
0.6667 
1.3333 
2.1000 

1.2222 
2.7778 
1.8889 
1.6000 
2.3333 
1.1000 
1.4444 
1.0 90Q 
2.0000 
1.4444 
2.3000 
1.00^0 
2.888^ 
2.6667 
2.2222 
2. 0^09 
1.1318 
1.5556 


0.4889 
0.7500 
1.1944 
0.6944 
0.9636 
2.2637 
1.7500 
0.2500 
0.5444 
G.9444 
3.1044 
0.3611 
1.3778 
0.7500 
1.4333 
0.7778 
G.C909 
0.2500 
0.7773 
3.7889 
0.5003 
8.3611 
2.2500 
1.4444 
0.6903 
0.5636 
0.2778 


ESTIMATED  PARAMETERS  FOR  NICKEL 


SERIAL  # 

516726 
515355 
516678 
515058 
509090 
524739 
516  3  0  9 
516094 
516949 
516750 
516575 
516534 
509081 
516592 
515931 
515603 
516380 
515515 
516731 
515447 
516048 
516063 
509064 

524646 
516627 
5161C1 
516132 


MEAN 

2.0000 
2.000^ 
3.3333 
1.4444 
1.0909 
2.9286 
1.0000 
2.3333 
3.2000 
2.5556 
2.8339 
3.4444 
0. 9000 
2.6667 
2.3000 
3.1111 
2.0000 
3.3333 
3.0000 
1.4000 
1.7778 
2.888° 
4.0000 
3.7778 
2.3636 
1.3091 
1.4444 


VARIANCE 

1.5  5  56 
2.00^1 
3.750") 
1.2778 
1.0909 
2.9945 
1.2500 
3.5000 
2.8444 
2.2778 
3.3611 
2.2773 
0.5444 
3.25O0 
1.5667 
3.6111 
3.4000 
4.0C00 
5.0000 
1.6000 
1.9444 
1.1111 
3.0000 
4.6944 
2.4545 
2.4903 
1.7773 
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REGRESSION  LIN<=S  FOR  ALUMINUM 

SEPIAL  M  INTERCEPT       SLOPE  RESIDUAL  SAMPLE 

VARIANCE    CORRELATION 

516726  8.7010       0.f>1562  2.5*15  ^.6->55 

515355  7.6667  0.13898  5.7870  0.5887 

516678  7.4444  0.11205  3.8638  0.2^41 

515058  7.444*  -0.11984  17.5090  -C.1703 

50909C  14.0000  0.13755  14.9953  0.5162 

524739  12.3571  0.00456  15.6695  0.0790 

5163C9  6.*444  0.07506  6.4149  0.7303 

5160^4  7.2222  0.13638  4.6821  0.4141 

l]^!  6.7000  -0.02389  4.8907  -0.3358 

516750  6.1111  0.00514  5.5361  0.0591 

516575  10.5556  0.03202  27.8192  0.1924 

516534  14.6667  0.02293  4.1051  0,3194 

5C9C81  7.0000  -0.01O15  25.4999  -0.0019 

5165^2  6.2222  0.02678  3.6684  0.3622 

515931  9.2000  0.04747  13.7188  0.4366 

5!56Ce  6.2222  -0.00061  7. 0731  -0.C07* 

516380  6.^-545  0.08528  16.5446  6.4035 

515515  13.8889  n. 08896  11.6034  °.5362 

516731  3.444*  -0.03564  11.9624  -0.3033 

515447  3.4011  0.11010  4.4045  0.74°5 

IJfS^i  §• 2^2°  0.0600*  4.1227  C.7311 

516C63  5.5556  0.03312  6.3846  C.3314 

5C9064  11.2222  0.02695  C.9913  0.9128 

524646  10.2222  0.05629  1.2945  0.8950 

516627  7.5455  0.06129  2.6692  0.7378 

516101  9.7273  0.02141  13.4872  0.3296 

516132  8.4444  0.04049  °.42O0  0.4659 


REGRESSION  LINES  FOR  IRON 

SERIAL  v  INTERCEPT       SLOPE  RESIDUAL  SAMPLE 

VARIANCE    CORRELATION 

516726  30.1000  0.06375  4.6265  0.9275 

515355  27.3333  0.11819  20.9135  ^.15^1 

|}!£Z§  Z$*155A  0.00537  27.4270  0l0348 

515058  22.6667  0.03220  30.0517  ol2093 

50909C  32.4545  C. 10673  2ol2778  ol327* 

524739  33.6^29  0.06821  26.7694  0.6716 

5163C9  25.4444  0.3*263  4C#2«*3  0^8695 

516094  26.6667  0.20229  615651  0.9034 

^169*9  24.7000  0.17268  3017814  ol7165 

516750  17.3333  0.13473  6*.6504  614137 

516575  35.5555  0.19532  16.51*1  ol8408 

516534  34.44*4  0.18827  Ulb030  018263 

5C9C81  25.5000  0.05879  3513315  0^359 

516592  17.4444  0.07297  812143  015779 

515931  30.1010  0.11599  13.5679  0  7662 

IJ5S28,  31.6667  3.02774  39.1967  0.1417 

51638C  24.3636  0.17684  65.1446  614236 

gl55X5  39.3333  0.16124  20.4983  0.6562 

5-6^31  21.1111  0.14852  17.8794  ol33*2 

515*47  29.2000  0.22241  2l!<U42  0^7159 

516C48  34.2222  0.2245*  92^5444  Q.b^ 

516C63  19.2222  ^.12472  30^4751  oltlAO 

IS22SJ  35. TTTS  0.07985  312291  0  9169 

52*6^6  40.1111  0.19645  13.1995  0.909Q 

Z\iHl  U'8132  0.18915  11.5735  oll^ol 

5161C1  29.8132  0.06576  69.5*54  0.*269 

516132  27.4444  0.10908  6.7836  o!8581 
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REGRESSION  LINES  FOR  COPPER 


SERIAL  « 


516726 
515355 
516678 
515C58 

524739 
516309 
516094 
516949 
516750 
516575 
516534 
5C9C81 
516592 
515931 
5156C8 
516380 
515515 
516731 
515447 
516048 
516C63 
509C64 
524646 
516627 
5161C1 
516132 


INTERCEPT 


14.5000 

11.5556 

11.7778 

7.5556 

9.6364 

10.1429 

8.8889 

9.2222 

9.6C0C 

10.7778 

10.0000 

11.4444 

5.7000 

8.8339 

10.1000 

15.7778 

9.*5*5 

11.2222 

10.3333 

10.2000 

9.7778 

10.1111 

14.0000 

11.6667 

10.6364 

15.6364 

3.4444 


SLOPE 


0.04810 
0.06015 
0.01226 
0.01855 
0.04157 
0. 03660 

0.07843 
0.02886 
0.06982 
•0.06335 
0.09021 
0.1^755 
0.01973 
0.063OQ 
0.0^496 
0.07357 
0.08577 
0.38609 
0.02090 
C. 06485 
0.05309 
C. 04695 
0. 0605 3 
0.07562 
0.09017 
0.04924 
0.03*03 


RESIDUAL 
VARIANCE 

4.6190 
5.8640 
4.043  5 
3.2891 
C.9376 
4.3039 
6.1235 

10.4862 
2.490  5 
5.5536 
2.9364 
6.O609 
4.66^5 
4.19S5 
1.7140 

11. 1454 

1C.0621 
1.3030 

11.0115 
3.2416 

10.5203 
6.8290 
1.6263 
2.0192 
1.1567 

25.3663 
2.4722 


SAMPLE 
CORRELATION 

0.8667 
0.7449 
0.2029 
0.3*93 
0.9364 
0.7714 
0.7527 
0.2314 
0.8250 
-0.5838 
0.8620 
0.7928 
0.5058 
0.6503 
0.7927 
0.5799 
0.4999 
0.8428 
0  .  1 9 1 C 
0.6137 
0.5102 
0.4339 
0.^639 
0.9074 
0.9255 
0.5052 
0.6538 
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